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^em  Marginal  equipment  dose,  rads (Si) 

Des  Equipment  susceptibility  threshold  dose,  rads(Si) 

Df(r,t)  Filter  dust  dose  rate  distribution,  rads(tissue)/hour-micron 

Df(r)  Filter  dust  dose  distribution,  rads(tissue)/micron 

D^(t)  Filter  dust  dose  rate  function,  rads(tissue)/hour 

Df(t)  Filter  dust  dose  function,  rads(tissue) 

DI(r,t)  Dust  inhalation  mass  rate  distribution,  grams/sec-micron 
Dl^(r)  Nasal  dust  inhalation  mass  distribution,  grams/micron 
Dlp(r)  Pulmonary  dust  inhalation  mass  distribution,  grams/micron 

Ojj  Dose  from  dust  in  nasal  cavity,  rads(tissue) 

Dp  Dose  from  dust  in  pulmonary  tract,  rads(tissue) 

Dpp  Dose  rate  at  an  electronics  piecepact,  rads(Si )/hour 

Djjj  Dose  rate  due  to  cockpit  suspended  dust,  rads(tissue)/hour 
D^(t)  Total  dose,  rads(tissue) 

Ejj  Relative  dose  effectiveness 

EDjjj  Effective  dose  for  10  percent  incidence  of  vomiting 

FE(r)  Filter  trapping  efficiency  function 

G Dust  mass,  grains 

ILBcm  Crew  member  initial  lung  burden,  uci/kg 

IV  Inhaled  volume  rate,  cm’/sec 

K,  Nondimensionalizing  constant,  gm(dust)/hour 

Kj  Nondimensionalizing  constant,  gm(dcst)/hour 

Kj  Nondimensionalizing  constant,  rads(tissue)/hour 
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K4 
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L 

Lf 

'-fo 

Mcp{'".t) 

«cp(t) 

%it) 

Mf(r.t) 

Mf(t) 

Mf(t) 

«t 

N 

NA(r) 

Pj(r) 

PT(r) 
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ABBREVIATIONS  AND  SYMBOLS  (cont'd) 

Nondimensional izing  constant,  rads(tissue)/’nour 

Nondimensionalizing  constant,  rads{tissua)/hour 

Nondimensional izing  constant,  gm(dust)/cm 

Cockpit  “"low  length,  cm 

Weapon  loading  factor,  megatons/kilometer^ 

Baseline  weapon  loading  factor,  megatons/kilometer^ 

Cockpit  dust  mass  rate  of  accumulation  distribution, 
grams (dust )/hour-mi cron 

Cockpit  dust  mass  distribution,  grams(dust)/micron 
Cockpit  dust  mass  distribution,  grams(dust)/micron 
Cockpit  dust  mass  function,  grams(dust) 

Engine  dust  mass  rate  of  flow,  grams(dust)/hour 
Engine  dust  mass,  gi  ,(ms(dust) 

Engine/filter  dust  mass  rate  of  flow,  grams(dust)/hour 
Engine/filter  dust  mass,  grams(dust) 

Filter  dust  mass  rate  of  flow  distribution,  grams/hour-micron 

Filter  dust  mass  distribution,  grams/micron 

Filter  dust  mass  rate  of  flow,  grams/hour 

Filter  dust  mass,  grams 

Total  mass,  grams 

Constant  of  proportionality 

Nasal  dust  accumulation  distribution 

Particle  settling  probability 

Pulmonary  dust  accumulation  distribution 

Cloud  radius,  kilometers 

Critical  settling  radius,  microns 
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ABBREVIATIONS  AND  SYMBOLS  (cont'd) 

Re  Reynol ds  number 

R(t)  Particle  size  fallout  (Nibbling  Mouse)  function,  microns 
TF  Cloud  exit  time,  hours  after  detonation 

TI  Cloud  entry  time,  hours  after  detonation 

Tmc  Time  of  mission  completion,  hours  after  detonation 
TSF  Time  scaling  factor 

T(t)  Shorthand  notation  for  a functional  evaluation  at  t 

TT(t)  Shorthand  notation  for  a functional  evaluation  at  t 

Horizontal  velocity,  cm/sec 
Vj.  Settling  velocity,  cm/sec 

Vj  Critical  settling  velocity,  cm/sec 

W Weapon  yield,  megatons 

a Length  measurement,  cm 

a^  Variable,  grams (dust)/cm^ 

b Constant 

bj  Cockpit  radius,  cm 

d Distance  measurement,  cm 

d,^  Distance  from  nasal  cavity  to  epigastric  region,  cm 

dj  Distance  from  pulmonary  tract  to  epigastric  region,  cm 

f(r)  Particle  number  distribution,  particles/cm’-micron 

f(r,t)  Particle  number  distribution  function,  particles/cm’-micron 

fj  Constant 

fj  Constant 

g Length  measurement,  cm 

h Cockpit  height,  cm 

1 Constant 
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ABBREVIATIONS  AND  SYMBOLS  (cont'd) 

m Constant 

'>’,0  Engine  air  mass  rate  of  flow,  gram/hour 

ac 

liiaf  Filter  air  mass  rate  of  flow,  gram/hour 
m^  i^^  mass  element 

m(r)  Mass  per  volume  distribution,  grams/cm^-micron 

n Constant 

p Radius  length,  cm 

r Particle  radius,  microns 

r^^  Particle  radius  defined  by  R(tp),  microns 

S^.  i^^  distance 

t Time,  hours  after  detonation 

t^  Cloud  exit  time,  hours  after  detonation 

t^  Cloud  entry  time,  hours  after  detonation 

t|_  Last  time  dunmy  variable  for  numerical  Integration 

t Particular  time  point  used  in  numerical  integration  to 

" generate  r^  via  r^^  = R(tn) 

X Cockpit  settling  distance,  cm 

y Plane  radius,  cm 

z Distance  above  plane,  cm 

iT  Cloud  penetration  time,  hours 

0 Angle  measurement 

6 Activity  breakpoint,  microns 

VI ' Mass  attenuation  factor 

Pj  Density  of  air,  grams/cm* 

Density  of  air  in  the  cockpit,  grams/cm^ 

dCp 

Pj(r,t)  Dust  density  distribution,  grams(dust)/cm^-micron 
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Pd(r) 

P.(t) 


ABBREVIATIONS  AND  SYMBOLS  (cont'd) 

Dust  density  distribution,  grams(dust)/cm^-micron 
Dust  density  function,  granis(diist)/cm^ 

Baseline  dust  density  function,  grams(dust)/cm^ 
Dust  particle  density,  grams(dust)/cm^ 

Cockpit  crossing  time,  sec 
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SECTION  I 
INTRODUCTION 

In  the  event  of  a nuclear  war,  strategic  aircraft  are  subject  to  accidental 
or  deliberate  penetration  of  the  late  time  (10  minutes  or  longer  after  detona- 
tion) radioactive  dust  clouds  generated  by  surface  detonations  of  large-yield 
nuclear  weapons.  Because  of  the  presence  of  numerous  clouds  from  enemy  first- 
strike  attacks,  the  majority  of  cloud  penetrations  will  occur  during  the  I -gh- 
altitude.  Continental  United  States  (CONUS)  e>;it  phase  of  the  strategic  mission. 

Penetrations  of  radioactive  clouds  by  manned  oircraft  pose  potential  hazards 
to  the  mission  of  the  aircraft  in  three  ways.  First,  the  crew  of  the  aircraft 
is  susceptible  to  ionizing  radiation  and  to  direct  contact  with  the  radioactive 
particles.  Second,  the  electronic  equipment  can  be  degraded  by  ionizing  radia- 
tion. Third,  the  dust  ingested  by  the  engines  could  cause  erosion  of  the  com- 
pressor blades  and/or  other  degradation. 

Consider  first  the  crew.  They  are  subject  to  ionizing  doses  from  immersion 
in  the  "infinite"*radiating  cloud,  from  dust  accumulations  on  the  aircraft 
exterior,  from  the  dust  accumulated  in  the  cockpit,  and  from  dust  accumulations 
elsewhere  in  the  aircraft.  Dust  accumulations  in  the  interior  of  the  aircraft 
are  due  to  the  ingestion  of  large  amounts  of  air  to  cool  heavy-duty  electronics 
equipment  and  to  provide  cockpit  pressurization  and  air  conditioning.  Because 
of  the  continuous  ingestion  of  outside,  contaminated  air,  the  crew  is  also 
subjected  to  direct  contact  with  the  radioactive  dust  and  to  inhalation  of  the 
dust  suspended  in  the  air.  All  of  the  factors  above  are  potential  crew-disabling 
hazards.  Each  must  be  thoroughly  investigated  to  determine  the  seriousness  of 
its  impact  on  crew  performance  and  to  determine  corrective  actions  required  to 
minimize  this  impact. 

Because  outside  air  is  used  to  cool  electronic  equipment,  radioactive  dust 
could  be  introduced  into  the  innermost  parts  of  the  electronics.  Although  the 
amount  of  dust  accumulated  in  any  one  Line  Replaceable  Unit  (LRU)  or  "black  box" 
is  relatively  small,  the  doses  accumulated  by  susceptible  piece,  parts  could  be 
quite  large  because  of  the  nearness  of  the  dust  to  piece  parts  requiring  cooling, 

♦"Infinite"  implies  a spherica’  radius  such  that  the  majority  of  photons  emitted 
outside  the  sphere  are  absorbed  b-  the  atmosphere,  i.e.,  about  3000  meters. 
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i.e.,  as  near  as  1 centimeter.  Therefore,  ionizing  dose  is  a potential  cause 
of  degradation  and  must  be  investigated.  (In  today’s  electronics,  the  piece 
parts  are  mounted  on  printed  circuit  cards  which  are  not  cooled  directly  by  the 
air.  The  air  flows  through  sealed  plenum  chambers  and  carries  away  heat  in  the 
wall,  generated  bv  high-power  devices.  There  is  no  direct  contact  between  the 
piece  parts  and  the  dust,  and  only  ionizing  gamma  radiation  is  of  concern.) 

The  potential  problems  caused  by  ingestion  of  the  dust  into  the  engine  will 
be  addressed  only  to  the  extent  that  the  amount  and  size  distrib.ition  of  the 
dust  will  be  determined.  This  information  could  be  used  by  engine  design  and 
test  engineers  uo  determine  whether  or  not  a problem  exists.  Detailed  investi- 
gation of  the  effects  of  the  Ingested  dust  on  engine  performance  is  not  within 
the  scope  of  this  report. 

The  first  stage  of  the  investigation  contained  in  this  report  is  the  deter- 
mination of  the  clrjd  characteristics  of  interest  and  the  formulation  of  simple 
engineering  model'  which  approximate  the  behavior  of  the  cloud  with  time.  Filter 
design  criteria  are  then  discussed.  (These  are  the  inputs  needed  by  a designer 
to  design  and  locate  the  filter  in  the  aircraft.)  Techniques  necessary  to  select 
a set  of  filter  design  criteria  are  presented.  Next  a detailed  examination  of 
the  cockpit  dust  accumulation,  the  effects  of  the  cockpit  dust  on  the  crew,  and 
the  application  of  the  filter  selection  techniques  to  choose  the  optimum  cockpit 
filter  are  presented,  A similar  examination  dealing  with  the  electronics  LRUs 
is  then  presented. 

In  the  detailed  analyses,  maximum  effort  was  expended  to  make  the  results 
general  and  applicable  to  any  mission,  threat,  aircraft,  and  electronics  equip- 
ment. The  results  are  generally  presented  in  nondimensional  form  where  the  non- 
dimensional  izing  constants  are  functions  of  the  mission,  threat,  etc.  However, 
the  calculations  of  the  dust  and  dose  accumulations  in  the  cockpit  depend  on 
specific  cockpit  size  and  geometry.  Therefore,  these  results  are  presented  for 
a four  to  six  man  cockpit,  representative  of  a manned  bomber,  transport,  tanker, 
and  other  large  aircraft.  If  cockpits  or  cabins  cannot  be  reasonably  represented 
by  this  model,  the  analysis  must  be  reaccomplished  with  a more  suitable  model. 

For  example,  the  model  and  results  for  the  cockpit  presented  herein  may  not  be 
adequately  representative  of  fighter  cockpits,  or  large  pressurized  cargo  areas 
of  transports. 
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SECTION  II 

SYNOPSIS  OF  CLOUD  CHARACTERISTICS 

1 . DUST  DENSITY  FOR  BASELINE  THREATS 

To  analyze  the  problems  associated  wif'  an  aircraft  penetrating  the  high- 
altitude,  dusty,  radioactive  clouds  associated  with  surface  detonations  of 
nuclear  weapons,  one  must  first  postulate  a reasonable  threat.  In  this  work, 
the  starting  point  is  the  results  obtained  by  Whitaxer  (ref.  1),  who  postula.jd 
two  separate  situations.  The  first  was  a single  surface  burst  of  a 4-megaton 
nuclear  weapon.  The  second  was  massive  simultaneous  surface  detonations  of  four 
hundred  5-megaton  weapons  distributed  uniformly  over  a 200  x 200-kilometer 
surface.  The  first  threai  may  be  likened  to  an  attack  on  a single  hard  target. 
The  second  may  be  likened  to  a first  strike  attack  on  USAF  missile  fields. 

Based  on  output  of  complex  dust  computer  codes  and  his  years  of  study  of 
nuclear  dust  clouds,  Whitaker  determineU  dust  densities  of  the  cloud  as  a func- 
tion of  time  and  the  spatial  extent  of  the  clouds  associated  with  the  above 
threats. 

Several  factors  should  be  understood  at  the  (-nset  which  bear  directly  on  the 
relation  of  the  dust  environment  to  specific  mission  and  threat  dependent 
variables.  First,  surface  bursts  maximize  the  dust  environments  and  for  this 
reason  they  are  emphasized  here.  Second,  attention  is  restricted  to  the  dust 
envi  onments  at  late  times,  i.e.,  aircraft  cloud  entry  10  minutes  or  more  after 
the  ourst.  This  restriction  is  primarily  based  on  the  practical  consideration 
that  at  early  times  the  aircraft  crew  would  be  more  concerned  with  the  immediate 
hazards  of  encountering  large  solid  particles  (centimeter  to  meter  range)  and  in 
interactions  with  severe  burst-induced  turbulence.  Third,  at  these  late  times 
the  major  part  of  the  burst-induced  turbulence  has  subsided,  very  large  particles 
have  fallen  out,  and  the  temperature  of  the  cloud  is  relatively  close  to  ambient. 
Therefore,  the  10-minute  point  marks  the  end  of  hydrodynamic  motion  and  the 
beginning  of  a phase  which  may  be  considered  pure  fallout  in  a comparatively 
static  atmosphere.  The  dust  cloud  is  generally  localized  around  the  burst  point 
and  has  a relatively  homogeneous  spatial  density  distribution  as  a result  of 
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prior  mixing.  Finally,  the  effects  of  any  spatial  density  gradients  which  exist 
within  the  cloud  are  assumed  to  be  of  little  importance  to  the  aircraft  because 
integral  effects  (e.g.,  total  mass  collected  and  associated  doses),  rather  than 
instantaneous  effects,  are  more  important  in  this  analysis. 

According  to  Whitaker,  at  10  minutes  after  detonation,  the  cloud  size*  has 
stabilized.  For  a zero  wind  situation,  the  cloud  size  based  on  Whitaker's  work 
and  Glasstone  (ref.  2)  can  be  approximated  by 

R = 5.7W®‘“  (1) 

where  R is  the  cloud  radius  for  t ^10  minutes  after  detonation  in  kilometers 
and  W is  the  yield  of  the  weapon  in  megatons. 

The  dust  density  obtained  by  Whitaker  for  the  single  4-megaton  detonation 
can  be  approximated  for  times  less  than  10  hours  after  detonation  by  the 
expression 


Pj(t)  » 4.33  X 10‘®(t-i*6  + 1.3U-0-5)  (2) 

where  o^  is  the  dust  density  in  grams  (dust)  per  and  t is  in  hours  after 
detonation. 

The  dust  density  corresponding  to  the  multiburst  case  can  be  approximated 
for  times  less  than  10  hours  after  detonation  by  the  expression 

0j|(t)  = 1.73  X 10-'^(t'‘-^  + 1.31t-°-3)  (3) 

It  is  noted  that  the  above  dust  density  functions  are  based  on  the  result  of 
detailed  calculations  by  comprehensive  computer  codes  which  include  fallout  as  a 
parameter. 

In  general,  the  dust  density  is  a function  of  the  number,  type,  yield, 
height  of  burst,  time  of  detonation,  spatial  distribution,  etc.,  of  the  weapons 
as  well  as  environmental  factors  such  as  the  composition  of  the  earth's  surface 
under  the  burst  points  and  the  atmospheric  conditions  in  the  vicinity  of  the 
burst  points.  In  this  analysis,  worst-case  factors  have  been  chosen  to  arrive 


*Cloud  is  assumed  cylindrical. 
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at  the  distributions  above.  These  assumptions  include  zero  wind,  representative 
soil  conditions,  dry  temperate  atmosphere,  and  1 megaton  of  dust  aloft  per 
megaton  yield  at  1C  minutes  after  detonation.  The  densities  above  were  obtained 
based  on  these  assumptions. 

2.  DUST  DENSITY  FOR  A GENERAL  THREAT 

The  dust  densities  ?bove  are  based  on  specific  threats.  To  support  a general 
analysis  addressing  other  threats,  a general  dust  density  is  required.  To 
realize  this  objective,  a weapon  loading  factor,  L^,  was  originated.  This  factor 
is  defined  as 


L ^ TOTAL  MEGATONS  OF  WEAPON (S) 
^ AREA  OF  DUST  CLOUD  IN  km^ 


(4) 


ibis  factor  is  an  indication  of  the  intensity  of  the  attack  and  is  related  to 
the  dust  density  because  of  the  assumption  that  1 megaton  of  dust  is  aloft  at 
10  minutes  after  detonation  per  megaton  of  yield. 

Consider  now  the  single,  4-megaton  weapon  detonation  as  a baseline  load 
factor,  L^  . 


nR2 


where 


(5) 


R =■•  5.7(4)'>-“ 


(6) 


Therefore, 


= 0.0125 


MT 

km^ 


(7) 


The  general  dust  density  can  be  expressed  in  the  form 

Lx 


(t) 


(8) 


or,  using  equations  (2)  and  (7),  this  expression  reduces  to 
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Pj(t)  = 3.46  X 10-6  L^(t-i*6  + 1.31t-°-3)  (9) 

As  a check  of  this  expression,  consider  the  multiburst  situation.  The 
weapon  loading  factor  is 


I 

L 


. 400  X 5 
f ■ 200  X 206 


= 0.05 


HT 

km^ 


The  dust  density  is 

Pd(t)  = 


or 

Pj(t)  = 


3.46  X 10-6  (o.05)(t-‘-6  + 1.3U-0-3) 
1.73  X 10-'^  (t-*-6  + 1.31  t-“"S) 


(10) 

(11) 

(12) 

(13) 


which  is  identical  to  the  result  previously  obtained  by  Whitaker. 

The  dust  density  for  any  given  threat  then  can  be  determined  by  first  calcu- 
lation of  the  weapon  loading  factor  and  substituting  into  the  general  expression. 
It  is  emphasized  again  that  this  entire  development  is  pertinent  only  for  surface 
detonations  and  megaton  size  weapons. 

3.  DUST  CLOUD  ACTIVITY 

Now  that  the  dust  density  is  determined,  the  specific  activity  must  be 
investigated.  It  is  assumed  that  fission  fragments  are  generated  at  a rate  of 
3 x 10^6  fission  fragments  per  megaton  (i.e.,  100  percent  fission  yield)  and  are 
uniformly  distributed  over  the  same  volume  as  the  dust  cloud  at  10  minutes  after 
detonation  (ref.  2).  Using  the  average  decay  relation  of  radioactive  debris 
from  Glasstone  (ref.  2)  and  Whitaker's  10-minute  result,  the  following  relation 
is  obtained. 


A(t)  = 4.22  X 10»2 


photons 

hour-gram(dust) 


where  t is  in  hours  after  detonation. 


04' 
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It  is  notad  that  since  the  yield  is  assumed  to  be  100  percent  fission, 
neutron  activation  of  surface  materials  is  not  a significant  factor  in  the 
cloud  ra^^ioactivity . If  a significant  portion  of  the  yield  is  fusion,  then 
fewer  fission  products  will  be  generated.  However,  because  of  the  higher 
energy  neutrons  released  from  the  fusion  process,  neutron  activation  of  surface 
materials  can  be  significant.  Therefore,  in  general,  the  cloud's  radioactivity 
is  not  a simple  function  of  the  yield  as  is  assumed  in  this  work,  but  is  quite 
complex.  However,  the  100  percent  fission  yield  assumption  should  yield 
results  which  provide  an  upper  bound  on  the  actual  results. 


I 
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SECTION  III 

COMPARISON  OF  DUST  CLOUD  MODEL  RESULTS 
I . DUST  CLOUD  MODELS 

The  expressions  for  the  dust  densities  and  specific  activities  developed  in 
the  previous  section  ./ere  used  In  the  development  of  three  cloud  models,  which 
approximate  the  actual  cloud  behavior.  These  cloud  models  were  used  In  a perfect 
filter*  analysis  of  an  aircraft  environmental  control  system.  In  this  section, 
the  various  cloud  models  developed  will  be  discussed  briefly,  and  the  results  of 
the  perfect  filter  analysis  for  each  model  will  be  compared.  Based  on  this  com- 
parison, one  dust  cloud  model  will  be  chosen  as  the  baseline  model  to  be  used  in 
the  subsequent  detailed  analyses. 

The  first  dust  cloud  model  developed  was  the  Uniform  Fallout  Model,  described 
in  detail  in  appendix  A.  This  model  Is  a first-order  approximation  to  the  actual 
dust  cloud.  It  takes  fallout  into  account,  but  In  a rather  simple  fashion.  The 
dust  cloud  in  this  model  may  be  considered  to  consist  of  dust  particles  uniform 
in  size  and  activity.  Therefore,  the  governing  equations  are  relatively  simple 
and  easy  to  use,  but  are  not  functions  of  particle  size.  An  obvious  disadvantage 
to  this  model  then  is  that  it  would  not  furnish  the  filter  designer  a trapped 
mass  distribution  as  a function  of  particle  size.  Because  of  the  lack  of  particle 
size  dependency,  the  results  may  also  not  be  realistic. 

The  second  model  developed  was  the  Improved  Fallout  Model  (appendix  B).  In 
this  model,  particle  size  was  taken  into  account  both  in  the  relation  for  the 
dust  density  and  the  specific  activity.  The  dens'‘ty  distribution  was  obtained 
through  the  consideration  of  a number-distribution  function  based  on  a typical 
soil  at  the  burst  point.  The  activity  as  a function  of  particle  size  was  assumed 
to  be  proportional  to  the  particle  volume  for  very  small  particles  and  to  the 


*A  perfect  filter  is  defined  to  be  a filter  which  traps,  or  collects,  all  dust 
particles  which  enter  the  environmental  control  system, regardless  of  the  size 
of  the  particles. 
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particle  surface  area  for  larger  particles.  This  assumption  is  based  on  the 
consideration  that  small  particles  are  predominately  condensed  bomb  debris, 
whereas  larger  particles  are  composed  of  an  inert  tr>re  of  soil  covered  with  a 
surface  coating  of  weapon  debris.  The  specific  activity  (particle  activity/ 
particle  mass)  is  a function  of  particle  size  and  time  after  detonation  (to 
account  for  the  radioactive  decay  of  the  fission  fragments).  Fallout  is 
modeled  after  Stokes  Law  which  states  that  for  small  particles,  i..“c  fallout 
velocity  is  proportional  to  the  square  of  the  particle  size.  This  model  is  a 
great  deal  more  representative  of  the  actual  cloud  behavior  but  is  still  a 
homogeneous  cloud  model,  i.e.,  no  variation  in  cloud  parameters  with  altitude 
I or  range  from  the  point  of  detonation. 

' The  third  and  last  model  developed  (appendix  C)  was  the  Zero  Fallout  Model. 

Tir'.s  model  represents  a limit  case  where  there  is  no  fallout  and  is  useful  in 
tins  report  only  as  an  extreme  worst-case  standard  of  comparison.  The  only 
* time  dependence  in  this  model  is  the  radioactive  decay  of  the  fission  products. 

The  dust  density  and  the  mass  rate  of  flow  to  the  filter  are  constant  during 
cloud  penetration.  The  use  of  this  crude  model  in  a penetration  analysis  would 
lead  to  gross  overdesign  of  the  filter  and  is  presented  only  as  a limit  case 
and  as  a standard  of  comparison. 

, It  is  noted  that  the  detailed  results  from  all  cloud  models  in  appendixes 

A through  C are  nondimensional . The  nondimensionalizing  constants,  Kj  through 
I . Kj,  (defined  in  equation  (A-44))  take  into  account  specific  threat,  aircraft, 

mission  profile,  distance  from  source  (when  applicable),  and  photon  energy. 

► ' It  is  also  noted  that  the  aircraft  is  assumed  to  enter  the  cloud  at  some 

cloud  entry  time,  TI,  and  exit  the  cloud  at  a time,  TF  .are  TF  = TI  + aT,  and 
aT  is  the  penetration  duration. 

2.  PERFECT  FILTER  RESULTS 

[ A major  objective  of  this  report  is  to  provide  results  which  can  be  used  to 

set  aircraft  filtration  requirements.  Because  of  the  necessity  to  provide 
adequate  protection  to  the  crew  and  electronic  equipment,  and  the  cost  and  weight 
restraints  on  any  aircraft  filter  design,  the  filter  must  provide  protection 
against  any  reasonable  worst-case  cloud  penetration,  but  must  not  be  grossly 
overdesigned.  The  cloud  model  then  must  be  capable  of  supporting  the  design 
criteria  selection  and  must  produce  realistic  results. 
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The  two  models  which  will  be  investigated  as  candidates  to  support  the 
penetration  analyses  are  the  Uniform  Fallout  Model  and  the  Imp»t)ved  Fallout 
Model,  It  is  apparent  at  first  glance  that  the  Zero  Fallout  Model  would  predict 
results  leading  to  gross  overdesign  of  the  f'^lter  and  is  not  a candidate. 

Consider  the  dust  mass  trapped  by  a perfect  filter  during  penetrations  of 
clouds  modeled  by  the  two  candidates.  The  total  mass  trapped  as  a function  of 
time  is  identical.  However,  the  Improved  Fallout  Model  has  the  capability  also 
of  predicting  the  particle  size  distribution  of  the  trapped  mass.  This  added 
capability  is  significant  and  would  provide  a designer  with  more  information  to 
aid  in  the  filter  design. 

The  most  significant  differences  in  the  two  candidate  models  are  in  the  dose 
and  dose  rate  predictions.  Figures  1 and  2 depict  the  cloud  immersion  dose 
rates  and  doses  respectively  as  a function  of  time  for  the  two  candidate  models 
with  the  Zero  Fallout  Model  results  shown  as  a standard.  The  Uniform  Fallout 
Model  dose  rate  decreases  with  time  at  early  times  at  a significantly  higher 
rate  than  the  Improved  Fallout  Model.  This  difference  is  caused  by  the  differ- 
ences in  the  cloud  models.  The  particle  sizes  and  associated  activities  in  the 
Uniform  Fallout  Model  are  all  equal.  Therefore,  fallout  equally  affects  the 
dust  density  and  the  cloud  activity.  The  dust  density  and  the  specific  activity 
of  the  dust  in  the  Improved  Fallout  Model  are  dependent  on  the  particle  size. 

The  larger  particles  are  less  radioactive  per  unit  mass  than  the  small  ones. 

Early  fallout  is  associated  with  the  large  particles  which  have  relatively  lower 
activities.  The  smaller,  more  radioactive  particles  remain  aloft.  Therefore, 
the  dose  rate  decrease  with  time  in  the  early  phases  of  the  cloud  is  less, 
although  the  overall  dust  density  behavior  with  time  is  identical.  At  later 
times,  the  Uniform  Fallout  Model  dose  rates  results  tend  to  level  out,  while  the 
Improved  Model  results  tend  to  decrease  at  a greater  rate.  This  difference 
again  is  attributable  to  the  same  model  parameter  differences.  In  the  Improved 
Fallout  Model,  the  smaller,  more  radioactive  particles  are  now  falling  out  of 
the  cloud  taking  with  them  relatively  greater  amounts  of  the  radioactive  debris. 

The  cloud  immersion  doses  for  TI  = 10  minutes  as  depicted  in  figure  2 also 
reflect  major  differences,  as  would  be  expected  since  the  dose  is  merely  an 
integral  over  time  of  the  dose  rate.  The  Improved  Fallout  Model  dose  results 
are  considerably  larger  in  value,  because  the  smaller,  more  radioactive  particles 
remain  aloft  and  cause  mere  dose  to  be  accumulated  for  a given  mass  of  dust. 
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The  doses,  for  TI  = 10  minutes  and  AT  = 120  minutes,  due  to  the  trapped 
filter  duse  are  depicted  in  figure  3.  In  this  case  also,  the  doses  predicted 
using  the  Improved  Fallout  Model  are  higher  than  those  predicted  by  the  Uniform 
Fallout  Model.  This  again  is  directly  attributable  to  the  presence  of  the 
smaller  particles  which  are  trapped. 

In  sumnary,  the  total  masses  collected  by  a perfect  filter  are  the  same  in 
both  models.  The  cloud  dust  density  as  a function  of  time  is  also  the  same. 

The  Improved  Fallout  Model  provides  additional  information  about  the  size  dis- 
tribution of  the  trapped  mass.  The  dose  rates  and  the  doses  predicted  by  the 
Improved  Fallout  Model  are  larger  than  those  predicted  by  the  Uniform  Fallout 
Model,  but  are  still  significantly  lower  than  the  limiting  results  from  the 
Zero  Fallout  Model. 

3.  SELECTION  OF  MODEL  FOR  ANALYSIS 

The  Improved  Fallout  Model  is  chosen  to  support  the  aircraft  penetration 
analyses  for  the  following  reasons: 

a.  More  detailed  results  may  be  obtained,  i.e.,  trapped  mass  and  filter 
dose  as  functions  of  particle  size. 

b.  The  model  more  closely  approximates  the  actual  cloud  behavior. 

c.  The  dose  results  predicted  by  the  use  of  the  Uniform  Fallout  Model  are 
not  conservative,  and  could  result  in  underdesign  of  the  filter,  which  could 
compromise  the  aircraft  mission. 

In  view  of  these  considerations,  all  further  analyses  will  be  based  on  the 
Improved  Fallout  Model  of  the  radioactive  dust  clouds. 
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Figure  3.  Model  Comparison,  Filter  Dust  Dose,  TI  * 10  Minutes 
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SECTION  IV 

FILTER  DESIGN  COK^1DERATIONS  i 

The  primary  purpose  of  filters  in  the  environmental  control  system  (ECS)  of 
an  aircraft  during  a general  war  mission  is  to  minimize  the  ionizing  dose  | 

accumulated  by  the  crew  and  electronic  equipment  during  radioactive  cloud  pene- 
trations. This  section  presents  (1)  the  filter  criteria  required  by  a designer  / 

to  design  a filter  and  position  it  in  the  aircraft;  (2)  the  procedures  necessary 
to  arrive  at  this  set  of  criteria  by  a planner;  and  (3)  the  procedure  for  adapt- 
ing the  perfect  filter  model  presented  in  appendix  B to  a real  filter.  In 
general,  the  protection  required  by  the  crew  differs  from  that  required  by  the 
electronic  equipment.  Therefore,  the  filter  design  criteria  would  differ. 

1.  DISCUSSIDN  OF  POINT  DESIGN  CRITERIA 

The  starting  point  of  a filter  design  criteria  selection  process  is  an 
operational  analysis  based  on  aircraft  basing  plans,  mission  routings,  mission 
profiles,  time  lines,  threats,  etc.  The  end  result  of  this  analysis  would  be 
the  determination  of  a Point  Design  Condition,  which  effectively  means  the 
fixing  of  the  parameters  TI,  TF,  L^,  and  aircraft  altitude.  (For  a given 
analysis,  the  aircraft  parameters  of  interest  would  also  be  known.)  One  major 
constraint  is  placed  on  this  result.  The  cloud  inmersion  ionizing  dose  accu- 
mulated for  the  Point  Design  Condition  must  not  exceed  the  susceptibility 
threshold  levels  of  the  crew  and  electronics.  Filters  cannot  alleviate  these 
doses,  and  if  they  are  exceeded, the  mission  is  in  jeopardy.  In  fact,  there 
must  be  some  margin  between  the  cloud  immersion  dose  and  the  susceptibility 
thresholds  to  give  some  flexibility  in  the  filter  criteria  choice.  For  example, 
if  no  margin  is  allowed,  the  filter  effectiveness  and  the  filter  location 
required  to  provide  the  necessary  protection  may  be  completely  unrealistic  and 
unattainable. 

If  the  Point  Design  Condition  is  not  realistic  in  view  of  the  above  limita- 
tion, then  an  iteration  is  necessary,  with  constraints  on  sone  of  the  operational 
variables,  so  that  an  acceptable  Point  Design  Condition  is  obtained.  This  report 
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will  not  dwell  any  further  upon  the  selection  of  the  Point  Design  Condition 
because  it  is  out  of  the  scope  of  this  work.  From  this  point  on,  it  is  assumed 
that  the  Point  Design  Condition  is  a known  input.  The  remaining  criteria  selec- 
tion process,  however,  will  be  examined  in  detail. 

2.  FILTER  DESIGN  CRITERIA 

To  design  a filter  and  locate  it  in  an  aircraft,  a designer  must  be  provided 
the  following  criteria: 

a.  Filter  trapping  efficiency. 

b.  Total  mass  of  dust  trapped  by  the  filter. 

c.  Mass  distribution  of  the  trapped  dust  as  a function  of  particle  size. 

d.  Minimum  separation  distances  from  the  filter  to  the  crew  and  electronic 
equipment. 

The  filter  trapping  efficiency  is  defined  as 

, _ dust  of  radius  r trapped  in  the  filter 
' total  dust  of  radius  r entering  the  til  ter  ' ' 

All  the  dust  which  is  not  trapped  by  the  filter  could  enter  the  cockpit  or  the 
plenum  chambers  of  the  electronic  equipment  and  become  a source  of  radioactivity 
which  could  adversely  affect  the  mission  completion  capability  of  the  aircraft. 

The  total  mass  of  dust  trapped  is  required  for  filter  sizing,  while  the  mass 
distribution  is  required  to  support  the  possible  necessity  to  stage  the  filter. 
For  example,  the  filter  may  be  a two-stage  design.  The  first  stage  may  trap 
large  particles  and  the  second  stage  may  trap  small  ones.  The  minimum  spacing 
criteria  would  provide  the  designer  with  the  input  necessary  to  develop  an 
envelope  of  satisfactory  filter  locations.  This  envelope  would  be  coordinated 
with  other  designers  to  determine  filter  locations  consistent  with  other 
functional  requirements. 

Selection  of  the  filter  criteria  is  based  exclusively  upon  the  marginal 
ionizing  dose  which  may  be  accumulated  by  the  crew  and  the  electronic  equipment 
due  to  all  sources  except  cloud  immersion.  The  differences  between  the  cloud 
immersion  dose  (which  in  the  proper  units,  i.e.,  rads(tissue)  or  rads(Si),  is 
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equally  applicable  to  both  the  crew  and  the  electronics)  and  the  susceptibility 
thresholds  of  the  crew  and  electronics  are  the  basis  for  the  filter  design 
criteria  selection. 

Based  on  available  information,  the  effective  dose  for  a 10-percent  incidence 
of  vomiting  (EDio)  in  man  is  about  70  rads(tissue)  mid-epigastric  dose  (ref.  3). 
This  dose  converts  to  about  100  rads(tissue)  incident  dose  which  is  appropriate 
to  the  work  at  hand.  Above  this  level,  nausea,  vomiting,  and  other  performance 
degrading  responses  appear  in  an  increasing  number  of  irradiated  humans.  The 
electronic  equipment  susceptibility  thresholds  are  equipment  dependent.  For 
existing  equipment,  the  values  are  the  thresholds  of  damage/upset  of  the  equip- 
ment to  ionizing  dose  and  are  obtained  through  an  assessment  of  the  equipment. 
Much  equipment  presently  under  development  has  a total  ioniring  dose  requirement 
laid  upon  it  by  the  responsible  Program  Office.  This  requirement  may  be  used  as 
the  susceptibility  threshold  in  this  case.  Otherwise,  more  detailed  analysis  of 
the  equipment  must  be  accomplished. 

At  this  point,  it  is  assumed  that  the  Point  Design  Condition,  which  yields 
the  cloud  itimersion  dose,  D^,  and  the  crew  and  equipment  susceptibility  threshold 
doses  (Ocs  and  Dgj)  are  known.  The  marginal  crew  dose,  and  the  marginal 
equipment  dose,  D^^,  are 


= 0 “0 
cs 

(16) 

s D - D 
es  c 

(17) 

These  doses  are  the  constraints  placed  on  the  crew  and  the  electronic  equip- 
ment filter  criteria  selection  process.  In  other  words,  the  total  dose  due  to 
the  trapped  dust  in  the  filters,  the  dust  accumulated  in  the  cockpit  and  avionics 
equipment,  and  the  dust  on  the  aircraft  exterior  must  not  exceed  the  marginal 
dose  levels. 

All  crew  doses  from  equipment  not  open-cycle,  air-cooled,  or  equipment  more 
than  several  meters  from  the  crew  may  be  neglected  in  the  crew  filter  criteria 
selection.  In  general,  for  today's  streamlined,  high-speed  aircraft,  the  amount 
of  dust  accumulated  on  the  aircraft  exterior  is  minimal  and  may  be  neglected. 
Because  of  the  generally  lower  vulnerability  of  the  equipment  than  the  crew,  the 
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only  significant,  dose  to  consider  in  the  equipment  filter  ci-itcria  selection  is 
that  due  to  dust  trapped  in  its  own  interior.  This  dust  would  accumulate 
primarily  in  the  cooling  plenum  chamber,  which  may  be  only  centimeters  away  from 
the  susceptible  devices  being  cooled.  Unless  the  equipment  is  located  immediately 
adjacent  to  a filter  with  a relatively  large  amount  of  trapped  dust,  the  plenum 
chamber  dose  will  far  exceed  all  other  sources.  For  crew  members,  the  only  sig- 
nificant dose  contributions  that  will  be  considered  here  are  the  doses  due  to 
cloud  imnersicn,  the  dust  trapped  in  the  filters,  and  dust  deposited  in  the  cock- 
pit. For  electronics,  the  only  significant  dose  contributions  are  due  to  dust 
deposited  in  equipment  which  are  cooled  by  an  open-cycle  coaling  system  and 
possibly  due  to  dust  deposited  in  nearby  filters. 

3.  REAL  FILTERS  VERSUS  PERFECT  FILTERS 


The  analysis  presented  in  appendix  B is  based  on  the  use  of  a perfect  filter 
in  the  environmental  control  system.  A perfect  filter  is  defined  as  one  which 
traps  all  dust  entering  the  filter  regardless  of  particle  size,  i.e.,  FE(r)  = 

1.0.  The  results  of  this  perfect  filter  analysis  were  used  in  the  previous  sec- 
tion to  compare  results  of  the  different  cloud  models.  As  will  be  shown  later, 
these  results  are  also  very  useful  in  a real  filter  analysis.  Of  course,  no  real 
filter  can  trap  all  input  dust  particles. 


Modification  of  the  perfect  filter  analysis  to  make  it  suitable  for  a real 
filter  is  straightforward.  Start  with  the  equations  for  the  mass  rate  of  flow 
of  dust  and  the  mass  of  dust  from  appendix  3 (i.e.,  equations  B-27,  B-28,  and 
B-29)  and  rewrite  them  to  include  the  filter  trapping  efficiency. 


Mf(r,t) 

Mf(r,t) 


Mf{r,t) 


= 9.9  X 10-2  FE(r)  r‘<>-5 

0.1  u <_r  <.  R{t) 

(18) 

= 1.9  X 10-2  r-0'5(t-ti) 

O.lu  <.  r i R(t) 

(19) 

= 9.9  X 10-2  FE(,  j r^-o-sCtn-ti) 

R{t)  5.  r^  s R(ti') 

(20) 

where  t^.  £ t^  <.  t,  and  r^  = R(tp). 
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From  this  point,  the  technique  used  to  determine  expressions  for  the  mass 
and  dose  accumulated  in  the  filter  is  identical  to  that  used  in  appendix  8 
(where  effectively  the  filter  efficiency  (FE)  was  1)  except  that  the  factor 
FE(r)  must  be  included  in  all  operations  upon  the  equations.  For  any  given 
F£{r)  the  same  technique,  which  consists  of  a mixture  of  analytical  and  numerical 
integration,  is  used  to  obtain  mass/dose  results  for  the  dust  trapped  in  the 
filter. 

For  simple  filter  efficiency  functions,  a simpler  graphical  technique  to 
obtai  this  result  is  appropriate.  This  technique  Involves  using  the  perfect 
filte  results  of  appendix  B.  For  example,  if 

FE(r)  = 0 O.lp  <.r  ;£.10u  (21) 

FE(r)  =1.0  r i.lOu 

use  the  perfect  filter  mass/dose  distributions  for  r 2.  lOy,  which  shows  the 
mass/dose  accumulated  for  a given  TI,  and  aT,  as  a function  of  particle  radius. 
The  total  mass/dose  is  obtained  from  the  cumulative  results  by  subtracting  the 
value  at  lOp  from  the  value  at  10‘*u.  More  detailed  examples  are  presented  in 
reference  4. 
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COCKPIT  ANALYSIS 

1.  POTENTIAL  CREW  HAZARDS  FROM  COCKPIT  DUST 

In  section  IV,  a general  procedure  which  led  to  the  determination  of  the 
filter  design  criteria  was  described.  One  of  the  key  elements  in  the  selection 
of  the  filter  criteria  was  the  dose  due  to  the  dust  not  trapped  by  the  environ- 
mental control  system  (ECS)  filter.  This  dust  enters  the  cockpit  and  is  a 
potential  hazard  to  the  crew.  This  section  presents  in  detail  the  determination 
of  the  ionizing  crew  dose  resulting  from  the  dust  entering  the  cockpit.  In 
addition,  the  effects  on  the  crew  due  to  inhalation  of  the  airborne  particles 
and  to  skin  contact  with  these  particles  are  investigated.  For  the  purposes  of 
illustration,  a Point  Design  Condition  and  a set  of  aircraft  parameters  will  be 
assumed  and  the  appropriate  filter  criteria  will  be  selected.  The  assumed 
conditions  were  chosen  deliberately  not  only  to  show  how  the  criteria  can  be 
obtained,  but  also  to  represent  a realistic  set  of  conditions  for  today's  cargo 
or  bomber  aircraft  so  that  this  example  may  be  directly  applicable  to  the 
majority  of  filter  criteria  selections  in  the  Air  Force  today. 

There  are  essentially  four  concerns  associated  with  the  dust  particles  which 
are  passed  by  the  filter.  The  first  is  the  consequence  of  the  airborne  dust  in 
the  cockpit.  It  is  a source  of  ionizing  dose.  The  second  concern  is  the  dust 
which  settles  out  in  the  cockpit.  It  is  also  a source  of  ionizing  dose.  The 
third  is  the  consequence  of  crew  inhalation  of  the  airborne  dust  particles.  The 
fourth  concern  is  the  direct  contact  of  the  radioactive  dust  particles  with  the 
skin  of  the  crew.  To  scope  the  problems  associated  with  cockpit  ingested  dust, 
zero  filtration  is  assumed  for  the  cockpit  ECS. 

2.  DUST  SUSPENDED  IN  COCKPIT  AIR 

The  dust  suspended  in  the  air  in  the  cockpit  will  be  examined  first.  This 
ionizing  dose  source  is  assumed  to  be  pertinent  only  during  the  actual  penetra- 
tion. The  air  flow  through  the  cockpit  is  assumed  to  be  sufficient  to  ensure 
that  the  dust  density  inside  is  the  same  as  that  outside,  corrected  by  the  cock- 
pit pressurization  factor,  It  is  also  assumed  that  the  dust  is 
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uniformly  suspended  and  that  the  cockpit  can  be  approximated  by  a sphere  of  1 
meter  radius.  The  dose  rate  at  the  center  of  the  sphere  then  due  to  the 
suspended  dust  is 


or 


1 00cm 


“"SO 


/ 


Pd(r,t)  A(r,t)  Cll  - FE(r)[  ^^e  ^Mpb'dbj 


4Tib2 


(23) 


^ ^cp  .,-lOOu'l 

Pg  L -I 


(24) 


Since  u'  is  approximately  lO'^  cm"',  the  factor  in  the  brackets  is  very 
small  and  results  in  a very  small  dose  rate.  Eecause  1-MeV  photons  are  assumed, 
the  constant  C = 4.88  x 10"'o  photon/cm^'  calculated  for  the  worst 

case,  i.e.,  TI  * 10  minutes,  aT  = 2 hours,  and  FE  = 0,  and  found  to  be  on  the 
order  of  1 rad(tissue).  Therefore,  this  source  of  ionizing  dose  is  negligible 
in  comparison  to  the  other  sources. 

3.  DUST  ACCUMULATED  IN  COCKPIT 

The  second  source  of  ionizing  radiation  is  due  to  the  dust  particles  which 
settle  out  in  the  cockpit  and  remain  in  the  cockpit  throughout  the  mission.  In 
addition,  those  particles  which  settle  on  the  crew  members'  clothing  and  nearby 
surfaces,  because  of  their  nearness  and  the  1/d^  nature  of  the  phenomenon,  can 
contribute  substantially  to  the  total  ionizing  dose  to  the  crew  members.  Thus, 
it  is  important  to  determine  the  amount  of  dust  settled  out  in  the  cockpit,  to 
determine  its  distribution  relative  to  the  crew  member,  and  to  determine  its 
contribution  to  the  total  ionizing  dose  accumulated  by  the  crew. 

To  obtain  meaningful  results,  the  actual  cockpit,  which  is  geometrically 
very  complex,  must  be  approximated  by  a modval  which  is  more  tractable  to  analysis. 
The  cockpit  is  modeled  by  a square  box  with  inlet  ducts  on  one  side  and  outlet 
ducts  on  the  opposite  side.  The  inlet  and  outlet  velocities  are  assumed  to  be 
uniform  and  on  the  order  of  the  uniform  velocity  of  air  across  the  box.  It  is 
also  assumed  that  no  settling  occurs  in  the  ducting  between  the  filter  and  the 
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inlet  duct  to  the  cockpit.  From  a mass  flow  continuity  equation  and  cockpit 
dimension  consideration,  we  can  calculate  a horizontal  velocity  of  the  air 
flowing  through  the  cockpit.  Neglecting  any  slippage  of  the  dust  with  respect 
to  the  air,  the  horizontal  velocity  of  the  dust  is  identical  to  that  of  the 
air  mass.  This  horizontal  velocity  of  the  dust  through  the  cockpit  can  be 
written  as 


where 


V = af 
''h  ■ iedO 

acp  cp 


(25) 


'’acp  ' density  inside  the  cockpit 

A,„  = cross-sectional  flow  area  of  the  cockpit 
cp 


Superimposed  on  the  dust  particles  is  the  previously  introduced  settling 
velocity,  Vj  (appendix  B,  equation  (B-2)),  which  is  a function  of  particle  size. 
The  time  required  for  a dust  particle  to  cross  the  cockpit  is 


(26) 


where  L is  the  length  of  the  cockpit  in  the  direction  of  the  flow.  During  this 
same  time  a particle  will  fall  a distance  of 


X 


= VsS  = 


(27) 


If  a particle  has  a settling  velocity  Vj  such  that  its  vertical  distance  of 
fall,  X,  is  equal  to  or  greater  than  the  cockpit  height  h,  then  it  will  be 
collected  in  the  cockpit.  Defining  a critical  settling  velocity,  V , as 

^s  = r''H 

any  particle  with  a settling  velocity,  V^,  greater  than  or  equal  to  Vj  will  be 
collected. 
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Since  it  was  assumed  that  the  initial  spatial  distribution  of  particles  is 
uniform  in  the  vertical  direction,  the  probability  of  a particle  with  < Vj 
settling  in  the  cockpit  will  be  a linear  function  of  its  settling  velocity.  For 
example,  if  a particle  has  Vj  = 1/2  7^,  it  will  only  fall  a distance  x = h/2  in 
the  time  it  takes  the  particles  to  move  through  the  cockpit.  Thus,  the  particles 
with  this  settling  velocity  which  are  initially  in  the  lower  half  of  the  cockpit 
will  be  collected  and  those  in  the  upper  half  will  not  h.ive  sufficient  residence 
time  to  reach  the  cockpit  floor. 

Defining  P5(r)  as  the  probability  of  a particle  settling  yields 


p5(r)  = 1 for  V5  > Vj 


(29) 


Ps(r) 


for  Vj  < 


(30) 


This  probability  can  be  related  to  particle  size  r in  microns,  because  Vj  is 
related  to  particle  size  through  Stokes  Law. 

Associated  with  the  critical  settling  velocity,  is  a critical  particle 
size,  R^,  given,  by 


and  the  related  probabilities  of  a particle  settling  in  the  cockpit  can  be 
given  by 


Pj(r)  = 1 .0  for  r >. 

86.4  r2 

Pc(r)  = -^P-'  ■ for  r < R 

' "af 


(32) 

(33) 


where  the  cockpit  volume  is  given  by  V^p  = L A^p.  It  is  noted  here  that  the 
actual  cockpit  volume  should  be  used  in  this  expression,  if  it  is  known.  The 
cockpit  height  can  be  chosen  as  the  maximum  height  if  more  conservative  results 
are  wanted,  or  the  effective  height,  which  is  the  actual  volume  divided  by  the 
horizontal  cross-section  area  of  the  cockpit. 
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It  is  noted  that  is  a "filter  efficiency"  of  the  cockpit.  Therefore, 

the  equations  for  the  mass  rate  and  mass  of  the  dust  accumulating  in  the 
cockpit  can  be  expressed  in  the  following  manner; 


M._(r,t) 

“ 9.9  X 10-2[l-FE(r)]  P5(r)r-o*5 

M.„(r,t) 

VP  - A n in-STi  n 


lu:S.r:iR(t)  (34) 


9.9  X 10-2[l-FE{r)]  Pj(r)r-®*5(t-ti)  ly  i r < R(t)  (35) 


= 9.9  X lO-Hl-FE(rn)]  Ps(rn)V“'^^’^n-ti)  R(t)  < ^R(ti)  (36) 


where  M (r,t)  is  the  rate  at  which  the  dust  settles  in  the  cockpit,  M (r,t)  is 
cp  cp 

the  mass  of  dust  accumulated  in  the  cockpit,  and  [1  - FE(r)]  accounts  for  the 
dust  collected  by  the  air  conditioning  in-line  filter.  A more  expanded  develop- 
ment of  the  equations  for  mass  deposition  in  the  cockpit  is  found  in  appendix  D. 

Calculation  of  the  dose  rates  and  doses  for  this  situation  is  not  as 
straightforward  as  the  calculation  associated  with  the  dust  trapped  in  the 
filter.  The  dust  settles  out  uniformly  over  the  entire  cockpit  and  is  not 
concentrated  at  a single  point.  Therefore,  to  obtain  dose  results,  some  sim- 
plifying assumptions  are  made.  Any  one  crew  member  would  be  subject  to  varying 
amounts  of  dose  from  each  dust  particle,  depending  upon  its  distance  from  the 
crew  member. 

From  a crew  dose  contribution  viewpoint,  the  effectiveness  of  an  element  of 
mass,  m.j , located  some  distance  s^  from  the  crew  member  is  proportional  to 
m.j/s.j^.  A relative  dose  effectiveness,  Ep,  is  defined  for  any  mass  distribution 


m. 
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If  the  entire  mass,  M^,  were  located  at  a distance  of  1 meter  from  the  crew 
member,  would  be  one.  Thus,  to  apply  the  results  (at  1 meter)  to  other 
situations,  multiply  them  by  the  relative  dose  effectiveness,  Ep,  determined  for 
the  dust  distribution  of  interest. 

In  this  work,  the  spatial  dust  distribution  assumes  that  98  percent  of  the 
dust  is  located  at  an  effective  distance  of  1 meter  (the  approximate  distance 
from  the  critical  mid-epigastric  region  of  the  seated  crew  member  to  the  cockpit 
floor  or  to  equipment  consoles)  and  that  2 percent  of  the  dust  is  located  at  an 
effective  distance  of  10  centimeters  (to  account  for  the  dust  deposited  on 
clothing  and  nearby  surfaces).*  These  percentages  were  estimated  from  horizontal 
surface  area  considerations  in  a typical  four-man  cockpit  and  result  in  a rela- 
tive dose  effectiveness  of  about  three.  This  value  has  been  used  in  the  cockpit 
dose  calculations  in  this  report.  (The  cockpit  dose  results  for  other  dust 
distribution  models  can  be  obtained  by  dividing  the  calculated  dose  (at  1 meter) 
by  three,  and  multiplying  by  the  relative  dose  effectiveness  for  the  dust  dis- 
tribution of  interest.)  A complete  mathematical  development  of  the  cockpit  dose 
is  presented  in  appendix  D.  Note  that  the  dose  ecuations  are  developed  based  on 
the  assumption  that  no  filter  at  all  is  installed.  Therefore,  all  the  dust 
enters  the  cockpit  and  FE(r)  = 0.  This  car*  '’s  called  the  imperfect  filter 
analysis  and  is  a limit  case,  just  as  the  perfect  filter  analysis  of  appendix 
B was  a limit  case  in  the  other  extreme.  It  will  be  shown  later  in  this 
section  that  these  two  limit  cases  can  be  used  to  great  advantage  in  the  filter 
criteria  selection  process  and  in  an  evaluation  of  existing  filters. 

For  the  purposes  of  illustration,  it  is  assumed  that  the  cockpit  effective 
height  is  2 meters,  the  cockpit  volume  is  1.83  x 10'^  cm^,  and  the  mass  rate  of 
flow  of  air  to  the  cockpit  is  1.24  x 10®  grams  per  hour.  These  parameters  should 
be  representative  of  the  cockpit  for  a crew  of  four  to  six  for  most  existing 
bomber  and  cargo  type  aircraft.  The  values  of  these  parameters  must  be  known 


‘This  rather  crude  distribution  was  used  because  it  is  impossible  a priori  with 
no  knowledge  of  the  specific  cockpit  geometry  and  equipment  location  to  determine 
the  exact  distribution.  This  distribution,  however,  should  be  a reasonable 
representation  (on  the  conservative  side,  deliberately).  The  equipment  in 
military  aircraft  is  arranged  on  consoles  within  easy  reach  of  the  crew.  Dust 
accumulations  on  these  surfaces  are  probably  1 meter  or  less  from  the  crew 
member. 
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or  assumed  in  order  to  calculate  the  critical  radius  of  the  dust  particles.  For 
these  values  R^.  = 20p.  The  cockpit  temperature  is  assumed  to  be  75°F,  the  cock- 
pit altitude  to  be  10,000  feet,  and  the  aircraft  altitude  to  be  30,000  feet. 

The  calculations  of  the  dust  mass  distribution  were  made  according  to  the 
equations  in  appendix  D,  and  the  results  are  presented  in  figures  01  through  09. 
Each  curve  provides,  for  a particular  cloud  entry  time,  the  dust  mass  distribu- 
tion function  as  a function  of  particle  size  for  various  penetration  durations, 
i.e.,  2,  5,  10,  15,  30,  60,  and  120  minutes.  The  entry  times  presented  range 
from  10  minutes  to  5 hours. 

The  results  show  that  for  r > R^  = 20y,  all  particles  are  collected  and,  in 
this  range  of  r,  the  results  are  identical  to  the  perfect  filter  results.  For 
r < 20y,  there  is  a marked  difference,  demonstrating  the  fact  that  small  particles 
tend  to  be  carried  through  and  out  of  the  cockpit  without  settling  out. 

The  next  series  of  figures  (appendix  D,  figures  DIO  through  D18)  displays  the 
cumulative  dust  mass  collected  in  the  cockpit.  The  cumulative  mass  is  defined  as 
the  contribution  to  the  cockpit  mass  attributable  to  all  particles  of  size  less 
than  or  equal  to  r.  These  figures  present  this  cumulative  mass  as  a function  of 

particle  size  for  a particular  cloud-entry  time  and  with  penetration  duration  as 

a parameter.  The  same  ranges  of  entry  times  and  penetration  durations  are  used. 
Observe  that  since  the  small  particles  tend  to  be  carried  out  of  the  cockpit 

with  the  air,  the  cumulative  mass  is  small  for  small  particles.  The  effect  of 

fallout  can  also  be  observed,  particularly  at  late  times,  because  the  lines  of 
constant  penetration  become  horizontal  in  the  larger  particle  size  range, 
indicating  that  no  particles  of  these  sizes  are  present. 

The  cumulative  mass  evaluated  at  r = 10,000  microns  is  equal  to  the  total 
mass  deposited  in  the  cockpit  for  any  given  TI  and  pT.  Figure  D19  depicts  this 
total  mass  accumulation  as  a function  of  time  after  cloud  entry.  A comparison 
of  figure  010  and  the  perfect  filter  results  for  the  same  TI  given  by  figure  833 
(adjusted  by  the  same  K^,  equation  59)  shows  that  an  appreciable  number  of  small 
particles  are  not  retained  in  the  cockpit.  This  difference  is  entirely  due  to 
the  fact  that  the  smaller  particles  are  ejected  from  the  cockpit  with  the  air. 

Since  the  smaller  particles  have  higher  specific  activities,  and  the  cockpit 
is  less  efficient  in  collecting  these  "hottest"  particles,  one  might  expect  that 
the  cockpit  dose  for  a given  TI  and  aT  would  be  less  than  the  perfect  filter 
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do5p  (at  1 meter).  However,  the  filter  dust  represents  a localized  source  of 
radioactivity,  whereas  the  settled  dust  in  the  cockpit  is  distributed  over 
horizontal  surfaces  within  the  cockpit  and  on  the  crew  members.  Because  of  the 
nearness  of  some  of  these  particles  to  the  crew  member,  the  particles  effective- 
ness from  a dose  standpoint  is  greatly  enhanced.  Therefore,  it  is  difficult  to 
state  a priori  which  dose  would  be  greater. 

The  cockpit  dose  results,  which  are  obtained  from  analytical/numerical 
integration  of  the  equations  developed  in  appendix  D,  are  presented  in  three 
series  of  graphs.  The  cloud  entry  times  and  the  penetration  durations  were 
deliberately  chosen  to  coincide  with  those  values  used  in  the  perfect  filter 
analysis  in  appendix  B.  The  first  series  of  figures  (figures  D20  through  D28) 
presents  the  cockpit  dust  dose  distribution  as  a function  of  particle  size  for 
a particular  cloud-entry  time  and  with  penetration  duration  as  a parameter. 

The  second  series  (figures  D29  through  D37)  presents  the  cumulative  cockpit 
dose,  i.e.,  the  cumulative  dose  due  to  all  particles  which  have  the  size  r or 
less,  as  a function  of  r for  a particular  cloud  entry  time  and  with  penetration 
duration  as  a parameter.  The  third  series  (figures  D38  through  D46)  presents 
the  total  cockpit  dose  due  to  all  the  dust  mass  which  has  settled  in  the  cockpit 
as  a function  of  penetration  duration  for  a particular  clouu  entry  time.  The 
pertinent  results  may  be  sealed  to  mission  lengths  other  than  the  one  presented 
by  techniques  identical  to  the  ones  discussed  in  appendix  B. 

It  is  apparent  from  these  results  that  the  ionizing  dose  from  dust  accumu- 
lated in  the  cockpit  for  the  no-filter  case  is  significant  and  must  be  accounted 
for  in  any  survivability  or  vulnerability  study. 

4.  INHALATION  OF  AIRBORNE  DUST  IN  COCKPIT 

The  above  discussion  investigated  the  ionizing  dose  which  the  crew  member 
would  accumulate  from  the  dust  ingested  into  the  cockpit.  This  dose  is  a result 
of  gamma  radiation  only.  The  next  concern  deals  with  the  inhalation  of  the 
airborne  dust  particles.  Host  modern  aircraft  have  a requirement  for  a "shirt- 
sleeve" cockpit  environment.  A shirt-sleeve  environment  means  that  the  crew 
member  would  be  comfortable  in  the  cockpit  with  normal  flight  clothing.  Tempera- 
ture and  pressure  are  maintained  such  that  no  special  equipment  or  clothing  is 
required  for  crew  comfort.  For  the  purposes  of  this  work,  it  is  assumed  that 
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no  gloves  are  worn  and  that  the  crew  member  is  breathing  cockpit  air.  Therefore, 
during  penetration  of  the  cloud,  the  crew  member  will  inhale  the  airborne  du;t 
and  retain  some  portion  of  it  in  his  body.  This  retained  dust  is  a source  of 
gamma  and  beta  radiation. 

The  total  dust  accumulated  in  the  respiratory  system  of  a crew  member  by 
breathing  air  which  is  contaminated  with  dust  will  now  be  estimated.  From 
reference  5,  the  fraction  of  particles  of  radius  r accumulated  in  the  nasal  and 
pulmonary  regions  may  be  approximated  by  the  following  relations. 

NA(r)  = -9.241  + 9.585  r»/2  - 52.52  r>/3  + 52.57  r'/"  (38) 

PT(r)  = 8.785e-'-‘*''  - 18.65e-2'”  + lO.ASe'^-^''  (39) 

The  above  functions  are  "filter  efficiencies"  of  the  nasal  cavity  and  the 
pulmonary  tract,  respectively.  For  a man  at  rest. 

Inhaled  volume/time  = (respiration/min)(tidal  volume) 

IV  = (15)(1000)  ^ 

IV  = 250  si 

The  total  mass  of  dust  inhaled  per  second  is 

0i(r,t)  = IV  pj(r.t)  ^ (43) 

Therefore, 

DI(r,t)  IV  p^j(r,t)  ^dt  (44) 

The  total  dust  trapped  in  the  nasal  and  the  pulmonary  tracts  then  as  a 
function  of  size  is 


(40) 

(41) 

(42) 
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Dl^(r)  = / IV  Pj(r.t)  -£E-NA{r)dt 

^ d 

Dlp(r)  = / IV  p^{r,t)  ^ PT(r)dt 
t".  d 


Nasal 


Pulmonary 


(45) 


(46) 


Here,  assume  that  the  crew  ionizing  dose  due  to  suspended  particles  in  the 
respi  atory  system  is  negligible  and  that  the  only  potential  dose  contributor 
is  the  dust  trapped  in  the  body.  The  “effective"  whole  body  dose  (due  to  gamma 
radiation)  associated  with  the  dust  trapped  in  the  nasal  and  pulmonary  tracts 
can  be  approximated  by  the  equations 


DIfj(n,t)  A(.-,t)dtdr 
4itd,.^ 


(47) 


•/-f 


A(r,t)dtdr 


4rd/ 


(48) 


where  and  Op  are  the  equivalent  doses  associated  with  the  nasal  cavity  and 
the  pulmonary  tract  and  d^  and  d^  are  effective  distances  from  these  dust 
accumulations  to  the  epigastric  region.  These  doses  are  determined  in  this 
manner  because  there  is  evidence  to  suggest  that  radiation  impinging  the  epi- 
gastric region  causes  nausea  and  other  symptoms  which  tend  to  cause  performance 
degradation  on  a short  term  basis  (ref.  3).  The  doses  calculated  previously 
were  incident  doses.  The  incident  dose  is  related  to  the  mid-body  and  whole- 
body  dose  by  the  relation 

Whole-body  dose  = 0.667  x Incident  dose  (49) 

Performing  the  indicated  operations  for  the  worst  case,  i.e.,  TI  = 10  minutes 
and  AT  - 120  minutes,  yields  a total  equivalent  incident  ionizing  dose  of  0.357 
rad(tissue).  This  dose  is  negligible  in  relation  to  the  other  doses. 
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The  inhaled  and  retained  dust  also  produces  physiological  damage  to  the 
tissue  of  the  respiratory  system.  This  damage  is  due  to  the  beta  particles 
associated  with  the  contaminated  dust  retained  in  the  body.  The  alpha  particle 
emission  associated  with  the  fission  products  rapidly  decreases  with  time  after 
detonation  so  that  even  at  10  minutes  after  detonation,  the  alpha  particle 
emission  is  negligible.  For  this  first  cut  approximation,  it  is  assumed  that 
the  gamma  ray  is  the  result  of  a.  disintegration  of  a fission  product.  This 
same  disintegration  produces  a beta  particle,  which  in  the  previous  discussions 
was  neglected  because  it  is  so  rapidly  absorbed  by  the  atmosphere.  For  this 
analysis,  however,  the  beta  is  quite  important  because  the  radioactive  material 
is  in  direct  contact  with  tissue  and  can  cause  damage  to  the  tissue.  The 
specific  activity  derived  in  appendix  B and  in  section  II  then  is  appropriate 
to  this  situation 


A.(t)  = (4.22  X 10>2)t-i’2  (50) 

P 

where  A^  is  in  beta  particles  per  hour-gram(dust).  In  terms  of  particle  size 


A,(r,t)  = 11.83  A,  0.1  i r < 20y  (51) 

P 1 P 

236.55  A,  t'»*2 

Ag(r,t) pi i i lOy  (52) 


where  Aj  is  defined  in  appendix  A. 


In  discussing  tissue  damage  due  to  beta  radiation,  different  units  are  used, 
i.e.,  the  microcurie.  Since  each  beta  particle  is  assumed  to  be  the  result  of 
a disintegration. 


® (3.7  X 10>»)(10-6) 


(53) 


with  the  unit  equation 
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disintegrations 

liour 

hour-gram(dust) 

seconds 

(54) 

di si nteqrat ions/second 

curie 

curie 

microcurie 

where  Cg  has  the  units  of  microcuries  per  gram  of 

dust,  and  3.7  x 10'“ 

disintegrations/second  is  the  standard  definition 

of  a curie.  Therefore, 

CJr,t)  = 3.75  X 10^  O.ly  < r < 20u 

P 

(55) 

20m  i r <.  lOS 

(56) 

The  total  number  of  microcuries  associated  with  the  dust  which  is  inhaled 
and  retained  by  a crew  member  can  now  be  determined. 

Although  there  is  a relatively  small  amount  of  information  available  in  the 
literature  about  humans  inhaling  radioactive  material  either  as  a result  of 
exposure  to  fallout  or  as  a result  of  incidents  involving  other  radioactive 
material,  the  conclusions  whicn  can  be  drawn  from  this  information  are  limited. 
The  reason  for  this  limitation  is  that  little  was  known  and/or  recorded  about 
the  precise  amount  inhaled,  the  activities  of  the  inhaled  material,  etc.. 
Therefore,  for  the  purposes  of  this  analysis,  the  human  data  are  of  little 
value  as  a standard  for  comparison. 

Because  of  the  impreciseness  associated  with  human  data,  other  areas  were 
explored  and  it  was  found  that  some  excellent  experimental  work  has  been  done 
by  the  Lovelace  Foundation,  Albuquerque,  New  Mexico  (ref.  5).  Carefully  con- 
trolled experiments  to  determine  the  effects  on  dogs  of  inhalation  of  various 
radioactive  aerosals  have  been  conducted.  All  of  the  experiments  were  tightly 
controlled.  The  dogs  were  observed  closely,  and  the  physiological  effects  were 
noted  by  professionals.  Therefore,  these  experiments  are  valuable  to  the 
present  work.  The  dog  used  was  the  Beagle,  which  is  smaller  than  a man  and,  of 
course,  is  not  identical.  However,  the  results  were  presented  in  terms  of 
Initial  Lung  Burdens  (ILB)  of  radioactive  material,  which  has  units  of  micro- 
curie/kg (body  weight).  The  body  weight  of  a warm-blooded  mammal  is  representa- 
tive of  its  vital  lung  capacity  because  the  lungs  furnish  oxygen  to  tissue,  and 
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the  more  tissue  to  be  supplied,  the  larger  the  lungs  required.  Therefore,  it  is 
assumed  that  the  Initial  Lung  Burden  which  can  be  calculated  for  a crew  member 
is  correlateable  (at  least  approximately)  to  the  Beagle  results  and  some 
conclusions  can  be  drawn  about  the  effect  of  the  inhaled  dust  on  the  human. 

The  Beagle  study  of  most  significance  here  is  the  one  (ref.  6)  wherein  a 
radioactive  isotope  in  fused  clay  was  inhaled.  The  isotope  has  a relatively 
short  half  life,  and  is  a beta  emitter.  These  properties,  as  well  as  the  fact 
that  fused  clay  was  used  as  a carrier,  are  similar  to  the  properties  of  the 
radioactive  particles  generated  by  surface  detonations  of  nuclear  weapons.  The 
Initial  Lung  Burdens  of  the  dogs  varied  from  80  to  5200  yCi/kg  and  the  associated 
lung  doses  from  1200  to  55,000  rads.  The  clinical  signs  observed  in  the  dogs 
that  died  were  those  typical  of  severe  pulmonary  injury  and  included  progressive 
increases  in  the  respiratory  rate,  abnormal  lung  sounds,  terminal  cyanosis  of 
mucous  membranes  and  increased  density  of  the  lungs  on  radiograph  examination. 

Although  the  nature  of  the  clinical  signs  did  not  differ  markedly,  the  time 
of  onset  of  the  symptoms  and  the  times  of  death  varied  considerably  with  the 
Initial  Lung  Burden.  The  dogs  with  high  ILBs  can  be  classified  into  three 
groups.  The  first  group  of  four  dogs  had  ILBs  of  2600  to  5200  uC1/kg  and  lung 
doses  of  32,000  to  55,000  rads(tissue).  Symptoms  were  first  observed  at  7 to 
10  days  after  exposure  with  deaths  at  7.5  to  47  days  after  exposure.  The  second 
group  of  18  dogs  had  ILBs  of  850  to  2400  yCi/kg  and  lung  doses  of  11,000  to 
29,000  rads(tissue).  Symptoms  were  first  observed  at  3 to  4 weeks  after 
exposure  with  a median  survival  time  of  85  days  after  exposure.  The  third  group 
of  17  dogs  had  ILBs  of  590  to  760  uCi/kg  and  lung  doses  of  7300  to  9400  rads. 
Initial  symptoms  were  observed  at  6 to  8 weeks  with  a median  survival  time  of 
185  days  after  exposure.  Dogs  with  ILBs  of  80  to  460  uCi/kg  and  lung  doses  of 
990  to  5700  rads(tissue)  showed  no  detectable  clinical  symptoms  up  to  735  days 
after  inhalation  exposure.  In  this  experiment,  no  potential  performance 
degradation  was  assessed.  Only  cliriical  observations  were  made. 

In  order  to  compare  ILBs  for  the  human  case  to  the  above  experiment,  one 
must  calculate  the  total  microcuries  trapped  in  the  lungs  of  a crew  member,  and 
divide  by  a body  weight.  The  total  amount  of  microcuries,  Cp,  associated  with 
the  dust  trapped  in  the  pulmonary  system  is 
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Cu  - 


(IV)  Pjj(r-.t)  ^PT(r)  Cg(r,t)dtdr 

d 


(57) 


Assume  that  an  average  crew  member's  body  weight  is  80  kg.  Therefore,  the 
crew  member's  Initial  Lung  Burden,  ILB^^,  in  pCi/kg,  is 


ILB 


_ Cu 
cm  SIS' 


(58) 


Performing  the  calculation  yields  (for  a TI  of  10  minutes,  a aT  of  120 
minutes,  the  same  mission,  and  aircraft  parameters  used  previously)  an  extremely 
worst  case  ILB  for  the  crew  member  of  325  pCi/kg.  This  level  is  relatively 
benign  when  compared  to  the  Beagle  levels  above.  Therefore,  because  this  level 
is  low,  it  can  be  predicted  with  reasonable  confidence  that  a crew  member  would 
experience  little  or  no  performance  degradation  over  the  duration  of  the  mission 
due  to  physiological  damage  caused  by  the  inhaled  dust.  (Although  long-term 
effects  may  be  significant,  here  only  the  time  associated  with  "bombs  on  target" 
is  pertinent.) 

As  noted  above,  the  actual  physical  damage  should  not  cause  performance 
degradation  on  a short-term  basis.  However,  the  psychological  effects  could  be 
significant.  An  individual's  performance  could  be  degraded  because  of  the 
emotional  reaction  to  the  awareness  of  the  potential  hazard  of  inhalation  of 
the  particles.  Therefore,  it  is  recommended  that  a warning  system  be  incor- 
porated in  tne  aircraft  to  notify  the  crew  that  the  aircraft  is  in  a radioactive 
cloud.  (Note  that  a warning  system,  and  not  a dosimeter,  is  recommended.  The 
dosimeter  readings  themselves  could  cause  psychological  problems.)  Upon 
initiation  of  the  warning,  the  crew  would  don  oxygen  masks  with  100  percent 
oxygen  for  the  duration  of  the  penetration.  This  simple  precautionary  action 
would  reduce  the  possibility  of  inhalation  of  the  particles  and  alleviate  the 
uneasiness  associated  with  the  inhaling  of  the  radioactive  particles.  This 
precautionary  measure  is  recommended  even  if  a cockpit  ECS  filter  is  installed 
because  no  filter  can  completely  remove  all  the  dust  particles  suspended  in  the 
ingested  air. 
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5.  EXPOSURE  OF  BARE  SKIN  TO  THE  RADIOACTIVE  DUST 

Because  of  the  large  number  of  variables  associated  with  this  potential 
hazard,  precise  calculations  of  the  s';in  dose  are  impossible.  For  example,  the 
crew  member  is  constantly  in  motion,  may  wipe  hands  on  clothing,  and  may  or  may 
not  be  perspiring.  All  these  factors  and  more  affect  the  results  of  a skin  dose 
calc  nation.  However,  an  upper  bound  for  the  dose  may  be  estimated.  This  upper 
bound  consists  of  estimating  the  skin  dose  by  assuming  all  of  the  dust  which 
would  settle  on  that  piece  of  skin  is  retained  on  the  skin  for  some  period  of 
time.  Bare  skin  is  also  assumed,  i.e.,  no  clothing  covering  the  skin  at  the 
point  of  exposure. 

Probably  the  main  region  of  concern  in  the  skin  would  be  the  dermis,  or 
corium.  (It  is  p"Otected  by  the  outer  layer  of  skin,  the  epidermis,  which  varies 
from  0.07  to  0.1?;  mm  over  most  of  the  body.)  The  corium  is  P to  4 mm  thick  and 
contains  numerous  nerve  endings  (ref.  7).  Beta-radiation  penetrating  the  corium 
could  cause  damage  which  would  be  sensed  by  the  nerve  endings  and  could  result  in 
severe  pain  to  the  crew  member.  Since  the  epidermis  is  about  0.1  mm  thick.,  it 
absorbs  all  beta  particles  with  energies  less  than  about  200  xeV  and,  therefore, 
provides  some  degree  of  protection  to  the  sensitive  corium.  Beta  particles 
whose  energies  range  between  200  keV  and  700  keV  are  absorbed  in  the  corium, 
and  beta  particles  with  energies  in  excess  of  700  keV  pass  through  the  cerium. 
(Assuming  that  the  skin  has  the  same  properties  as  water,  the  graphs  in  refer- 
ence 8 were  used  as  the  basis  for  these  statements.)  As  a worst  case,  it  is 
assumed  that  all  the  beta  particles  emitted  by  the  radioactive  dust  which  settle 
onto  the  skin  possess  energies  in  excess  of  700  keV.  Although  this  assumption 
is  somewhat  conservative,  study  of  the  energi  s of  beta  particles  emitted  by 
typical  deca,  'ng  fission  products  Indicate  that  the  majority  of  the  beta 
particles  art  quite  energetic.  Therefore,  this  assumption  is  not  totally 
unrealistic. 

From  reference  7 it  is  noted  that  the  energy  loss  of  a beta  particle 
impinging  on  any  substance  is  approximately  2 MeV/gm/cm‘.  Therefore,  the 
conversion  constant  from  beta  particles  per  centimeter^-  to  rads  (tissue), 

Cj,  is 
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S 


100  / rads(tissue) 


C , 32  X lO-'*  rads  (tissue) 
6/cm^ 


(59) 


If  it  is  now  assumed  that  each  photon  generated  by  a disintegration  is 
accompanied  by  a beta  particle,  then 


A(r,t)  = Ag(r.t)  (60) 

lere  A.  is  the  specific  activity  of  the  radioactive  dust  in  8/gm(dust)-hour  and 

P 

A(r,t)  is  defined  in  section  II.  It  is  now  assumed  that  the  dust  settles  out  on 
the  exposed  skin  in  a uniform  manner  and  that  half  of  the  beta  particles 
released  by  the  layer  of  disintegrating  fission  products  penetrate  the  skin  in 
a normal  direction.  The  other  half  escapes  in  the  opposite  direction  and  is  not 
of  concern.  The  approximate  skin  dose  then  is 


Ds(r.t) 


(61) 


Cg  Hcp(r.t)  Ag(r.t) 


dt 


(62) 


It  is  noted  that  this  expression  differs  only  by  a constant  from  the  dose 
from  dust  accumulated  in  the  cockpit,  i.e.,  equation  (21).  In  other  words 


Ds(r.t)  D^.p(r,t) 
•'a  K, 


(63) 


where 


C Aj  r.i^,  a.j 
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l ' = ‘'S  ”*af 

' ' Oa<2H3)  a;p 

where  the  factor  of  2 in  the  denominator  of  Kj  accounts  for  the  assumption  that 
only  half  of  the  beta  particles  penetrate  the  skin,  the  factor  of  3 is  associated 
with  the  assumed  cockpit  dust  distribution  and  is  not  appropriate  for  the  skin 
dose  calculation.  The  factor  A^p  is  the  area  of  the  horizontal  portion  of  the 
cockpit,  and  a.  and  Aj  are  constants. 

Performing  the  indicated  operations  yields  the  relation 

D5(r,t)  = 15.0  Dpp(r,t)  (64) 

Therefore,  the  skin  dose  results  may  be  obtained  from  figures  D-20  through  D-47 
simply  by  multiplying  these  results  ty  15.0.  Note  that  the  skin  dose  results 
are  quite  significant  even  for  rather  short  penetration  durations. 

The  severity,  time  of  appearance,  and  duration  of  skin  response  as  a function 
of  radiation  depends  on  a complex  interaction  of  many  factors,  e.g.,  dose  rate, 
anatomical  \egion  exposed,  area  of  skin  irradiated,  presence  of  other  irritants 
(e.g.,  clothing  chafing  the  irradiated  area),  energy  of  impinging  beta  particles, 
and  depth  of  penetration  into  tlie  skin  {see  reference  3 for  a more  detailed 
discussion).  However,  for  the  case  assumed  above,  i.e.,  beta  particles  with 
energies  in  excess  of  700  keV,  the  surface  dose  required  to  produce  recognizable 
transopidermal  injury  is  about  2000  rads (tissue).*  Therefore,  unprotected  skin 
in  a cockpit  which  is  supplied  with  unfiltered  air  during  cloud  penetrations  of 
relatively  short  durations  could  be  a performance  crippling  factor.  Some  degree 
of  protection  would  be  afforded  the  crew  member  by  wearing  gloves  and  a helmet 
during  the  mission.  However,  reliance  on  protective  clothing  alone  would  violate 
the  cockpit  "shirt-sleeve"  environment  requirement. 


‘This  is  the  dose  required  for  50  percent  of  the  human  population  to  experiepce 
visible  skin  damage.  The  human  susceptibility  threshold  is  probably  lower, 
i.e.,  on  the  order  of  800  to  1000  rads  (tissue). 
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6.  EXAMPLE  COCKPIT  FILTER  SELECTION 

In  the  filter  criteria  selection  process,  the  perfect  filter  results  of 
appendix  B are  used.  To  use  these  results,  the  nondimensionalizing  constants, 
Kj  through  must  be  determined.  The  altitude,  the  aircraft  characteristics, 
i.e.,  cockpit  size,  mass  flow  rates  of  air  to  the  cockpit,  etc.,  have  been 
assumed  earlier  in  this  section.  Now  assume  that  the  threat  corresponds  to  the 
massive  multiburst  situation,  so  that  L^  is  fixed. 

The  pertinent  constants  (equations  (A-44))  then  are  determined  to  be 


= 480.76  9Bi^ 

(65) 

K3  = 7.878 

(66) 

= 12.724  radsl^ssuel 

(67) 

In  the  following  discussion,  the  values  obtained  for  the  perfect  filter 
cases  arc  the  values  read  off  the  ordinate  of  the  particular  graph  multiplied 
by  the  pertinent  constant. 

Now  that  a cockpit  deposition  model  has  been  developed,  and  the  primary 
short-term  concern  associated  with  the  dust  which  enters  the  cockpit  has  been 
determined  to  be  the  ionizing  crew  dose  from  the  nearby  dust  accumulations, 
filter  design  criteria  can  be  chosen.  This  selection  procedure  is  based  solely 
upon  the  necessity  to  protect  the  crew  from  excessive  accumulation  of  ionizing 
dose. 

It  is  assumed  for  the  purpose  of  this  exercise  that  the  Point  Design  Condi- 
tion has  been  fixed  and  leads  to  the  case  of  TI  = 30  minutes  and  aT  = 10  minutes. 
From  the  cloud  immersion  dose  results  {figure  B-61),  the  cloud  immersion  dose  is 
47  rads(tissue).  Since  it  was  stated  earlier  that  the  EDj^  for  vomiting  was 
100  rads(tissue)  incident,  there  is  a margin  of  53  rads.  In  other  words,  the 
total  crew  dose  from  the  cockpit  dust  and  the  filter  dust  must  not  exceed  53 
rads.  It  is  assumed  that  all  other  sources  of  ionizing  dose  due  to  cloud 
penetration  are  vanishingly  small. 
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To  alleviate  all  potential  hazards  from  the  cockpit  dust,  23  rads(tissue) 
are  chosen  as  the  marginal  cockpit  dust  dose,  and  30  rads(tissue)  as  the  mar- 
ginal filter  dose.  Going  to  figure  031  and  choosing  aT  = 10  minutes,  it  is 
found  that  the  23  rads(tissue)  limit  corresponds  to  a particle  size  of  18 
microns.  Therefore,  the  filter  trapping  efficiency  must  be 


FE(r)  = 1.0 

r >.  18u 

(68) 

FE(r)  = 0.0 

r < 18p 

(69) 

Now  go  to  figure  835,  which  is  the  perfect  filter  dose  (at  1 meter).  The 
dose  for  TI  - 30  minutes  and  aT  = 10  minutes  is  71  rads(tissue).  However,  the 
particles  smaller  than  18  microns  are  passed.  The  dose  associated  with  particles 
18  micron  and  smaller  is  35  rads.  Therefore,  the  total  filter  dose  is  71  - 35  = 
36  rads(tissue)  at  1 meter.  Therefore,  to  attenuate  this  dose  to  the  marginal 
30  rads(tissue),  the  filter  must  be  located  at  a distance  greater  than  1 meter 
from  the  crew  member.  The  distance,  d,  corresponding  to  the  30  rads(tissue)  is 

^ = 30  (70) 

d2 

d = 1.2  m (71) 

Therefore,  the  filter  must  be  located  at  least  1.2  meters  away  from  the  nearest 
crew  member. 

The  only  remaining  part  of  the  design  criteria  needed  is  the  total  mass, 
and  its  distribution  with  size,  trapped  by  the  filter.  To  obtain  this  informa- 
tion go  to  figures  87  and  B16.  From  figure  87  the  mass  distribution  is  obtained 
by  recalling  that  all  particles  less  than  18  microns  are  passed.  The  remainder 
of  the  particles  are  trapped.  Therefore,  the  mass  distribution  is  represented 
by  the  aT  = 10  minute  curve  with  the  distribution  function  below  18  microns  set 
identically  to  zero.  The  total  mass  trapped  by  the  filter  can  be  obtained  from 
figure  816.  Read  the  mass  collected  for  the  aT  = 10  minute  curve  corresponding 
to  r = 10,000g,  which  is  326  grams.  Subtract  from  this  value  the  mass  corre- 
sponding to  r = 18y,  which  is  72  grams.  The  result,  254  gm,  is  the  total  mass 
accumulated  by  the  filter. 
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Therefore,  it  is  apparent  that  given  the  Point  Design  Condition,  which  is 

4 

based  on  an  operational  analysis,  it  is  quite  straightforward  to  choose  filter 
design  criteria  which  provides  the  crew  necessary  protection  from  the  radio-  ^ 

active  dust.  This  design  criteria  may  be  specified  to  a designer,  who  can  design  ■! 

the  filter  and  locate  it  in  the  aircraft  without  any  knowledge  of  nuclear  surviv-  i 

ability/vulnerability  (S/V)  or  nuclear  environments. 

This  method  may  also  be  used  to  evaluate  existing  filters.  The  only  differ-  ‘ 

ence  would  be  that  the  trapping  efficiency  is  now  fixed  and  the  associated 

cockpit  dose  is  not  variable,  but  also  fixed.  The  filter  dose  then  fixes  the  / 

minimum  distance  to  the  crew.  The  filter  mass  as  dictated  by  the  filter 

efficiency  can  be  determined  and  compared  to  the  capacity  of  the  existing  filter.  ' 

I ' 

i 


► 
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SECTION  VI 

ELECTRONIC  EQUIPMENT  FILTER  ANALYSIS 

1.  GENERAL  CONSIDERATION 

In  the  previous  section,  a detailed  analysis  was  presented  which  resulted 
in  design  criteria  for  the  cockpit  environmental  control  system  filter.  Earlier 
it  was  stated  that  separate  filters  would  probably  be  required  for  the  cockpit 
and  the  electronic  equipment  because  of  different  filtration  requirements. 
Therefore,  it  is  necessary  now  to  investigate  the  electronic  equipment  protec- 
tion requirements  and  to  select  an  appropriate  set  of  filter  design  criteria. 

In  one  respect,  the  electronic  equipment  analysis  is  simpler  than  the  cock- 
pit analysis.  There  are  no  biological  considerations  other  than  ensuring  that 
the  crew  does  not  accumulate  excessive  dose  from  the  dust  trapped  in  each  line 
replaceable  unit  (LRU)(a  self-contained  "black  box").  Generally,  the  amount  of 
dust  trapped  is  relatively  small  and  the  LRUs  are  located  in  equipment  bays 
considerable  distances  away  from  the  aircrew.  Therefore,  the  effects  on  the 
crew  generally  are  negligible.  However,  in  another  respect,  the  analysis  is 
more  complex.  Each  LRU  in  the  aircraft  has  a potentially  different  power 
requirement  and  dissipation.  Therefore,  each  could  require  a different  amount 
of  cooling  air  and  each  could  have  a different  plenum  chamber  geometry.  There- 
fore, the  analysis  must  consider  each  LRU  Independently.  The  filter  design 
criteria  then  would  be  that  necessary  to  provide  the  protection  required  by  the 
most  susceptable  LRU. 

2.  PRELIMINARY  ANALYTICAL  CONSIDERATIONS 

Since  the  number  of  LRUs  in  a modern  aircraft  is  quite  large,  some  initial 
considerations  to  reduce  the  number  to  be  analyzed  are  in  order.  First,  only 
those  LRUs  which  are  mission  critical  need  be  considered.  Line  replaceable 
units  are  mission  critical  if  their  operations  are  necessary  to  the  successful 
completion  of  the  mission.  Second,  only  those  LRUs  which  are  air  cooled  by 
open  cycle  flow  techniques  must  be  considered.  Much  of  the  equipment,*  such  as 


*0nly  cloud  immersion  doses  are  pertinent  to  these  LRUs. 
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controls  and  displays,  and  other  low  heat  dissipation  LRUs,  do  not  require 
special  cooling.  Other  LRUs  are  cooled  by  a closed  cycle,  or  a regenerative 
cooling  technique.  This  type  of  cooling  technique  reuses  the  same  working 
fluid,  air.  Freon,  etc.  This  fluid  is  forced  through  the  LRUs,  through  a heat 
exchanger  to  remove  excess  heat,  and  then  back  through  the  LRUs.  Therefore, 
little  or  no  outside  air  is  introduced  into  the  LRUs.  The  open  cycle  cooling 
technique  uses  engine  bleed  air,  cools  it,  forces  it  through  the  LRUs,  and  then 
ejects  it  overboard.  Therefore,  there  is  a continuous  flow  of  outside  air 
through  the  LRUs  which  results  in  dust  deposition  in  the  plenum  chamber.  Third, 
only  those  LRUs  using  components  which  are  susceptible  to  total  ionizing  doses 
need  be  considered.  For  example,  those  LRUs  based  on  vacuum  tube  technology 
will  not  be  affected  by  total  dose  accumulations.  Only  the  latest  state-of-the- 
art  semiconductor  components  such  as  metal  oxide  semiconductors  (MOS)  devices 
are  susceptible. 

3.  LINE  REPLACEABLE  UNIT  ANALYSIS 

After  the  list  of  LRUs  requiring  analysis  has  been  reduced  to  a minimum,  the 
remaining  LRUs  must  be  analyzed.  The  first  step  in  this  analysis  is  the  develop- 
ment of  a dust  deposition  model  for  the  LRU  plenum  chamber.  The  basic  mechanism 
of  deposition  is  the  same  for  the  LRU  as  it  was  for  the  cockpit.  Therefore,  the 
same  model  will  be  used.  The  critical  particle  size,  R^  (defined  in  equation 
(31))  Isof  particular  Importance.  This  parameter  is  a function  of  the  plenum 
chamber  geometry  and  the  mass  rate  of  flow  of  the  cooling  air.  Both  these 
variables  could  vary  from  LRU  to  LRU  and  must  be  considered  in  each  LRU  analysis. 

The  approach  taken  at  this  point  is  to  obtain  general  results  which  can  be 
applied  to  any  specific  LRU.  Recall  that  in  the  cockpit  analysis  a representa- 
tive example  was  presented  as  an  illustration.  The  critical  particle  size  was 
20  microns.  The  LRU  critical  particle  size,  however,  may  vary  considerably  from 
LRU  to  LRU.  Therefore,  the  LRU  critical  particle  size  is  an  Important  parameter 
in  the  general  results. 

To  be  precise,  the  same  type  graphical  results  presented  in  the  cockpit 
deposition  analysis  should  be  presented  here  for  each  and  every  critical  particle 
size.  However,  presentation  of  these  results  would  require  thousands  of  graphs. 
Therefore,  another  approach  will  be  investigated.  In  the  cockpit  analysis,  the 
margin  between  the  crew  susceptibility  threshold  and  the  immersion  dose  was  tens 
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of  rads  (tissue).  Here  the  margin  between  the  LRU  susceptibility  and  the  iiiinrr- 
sion  dose  is  generally  thousands  of  rads(Si).  (It  should  be  noted  that  rads(Si) 

= 0.922  rad(tissue)  for  the  assumed  1 meV  photons.)  In  this  analysis,  extreme 
precision  and  accuracy  are  not  so  critical.  Therefore,  a simplified  approach 
will  be  described  which  yields  general,  slightly  conservative  results,  and  which 
can  be  simply  applied  to  any  LRU. 

The  equations  which  were  developed  in  appendix  D for  the  cockpit  are  appro- 
priate to  this  analysis  if  the  cockpit  critical  particle  size  of  20  microns  is 
replaced  by  the  general  parameter  R^  and  all  of  the  constants  are  divided  by  3 
to  account  for  the  magnification  factor,  E^,  used  earlier  in  the  cockpit  analysis. 
For  any  given  R^  these  equations  were  analytically/numerically  integrated  using 
the  Air  Force  Weapons  Laboratory  (AFWL)  CDC  6600.  The  results  of  interest  in 
this  simplified  approach  are  the  total  mass  of  dust  collected  in  the  LRU  and  the 
total  dose,  at  1 meter,  associated  with  the  total  mass  which  accumulates  over 
the  mission.  In  other  words,  in  this  analysis  the  distributions  of  the  mass  and 
dose  with  respect  to  particle  size  are  not  of  interest.  The  total  mass  collected 
can  be  used  to  determine  the  volume  taken  up  by  the  dust  in  the  plenum  chamber 
to  support  plenum  chamber  clogging  investigations.  The  total  dose  will  be  used 
to  establish  the  total  dose  incident  on  susceptible  piece  parts  in  the  LRU. 

It  is  apparent  that  the  Integration  of  the  equations  for  each  R^,  yields  a 
single  data  point.  Numerous  values  of  must  be  chosen  so  that  the  results  are 
comprehensive  enough  to  support  general  LRU  analysis.  These  operations  have 
been  performed  and  the  results  are  presented  in  figures  4 through  21.  Each 
figure  represents  a single  different  cloud  entry  time  and  each  curve  on  each 
figure  represents  a single  different  penetration  duration.  The  ordinate  is  the 
nondimensional  mass  or  dose  used  in  appendix  B and  the  abscissa  is  the  critical 
particle  size.  The  dose  curves  are  30-hours-after-detonation  cases.  The  doses 
for  other  times  can  be  obtained  by  use  of  the  time  scaling  factor  discussed  in 
appendix  B.  It  is  noted  that  each  point  on  these  curves  corresponds  to  the  area 
under  mass  and  dose  distributions  similar  to  those  presented  in  the  cockpit 
analysis  section. 

The  doses,  at-l-meter,  are  now  available.  However,  these  doses  were  deter- 
mined based  on  an  assumed  point  source  of  radioactivity.  The  dust  which  settles 
out  in  a plenum  chamber  is  not  a point  source.  It  is  distributed  as  a planar 
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source.  Therefore,  the  point  source  results  must  be  modified  to  take  this 
difference  into  account.  A plenum  chamber  model  will  be  introduced  which 
simplifies  the  complex  plenum  chamber  geometry,  and  from  which  a correction 
factor  can  be  devised.  It  is  assumed  that  the  dust  trapped  is  uniformly 
Hi'stri hiited  nupr  a n'rrnlar  area,  equal  to  the  plenum  chamber  area.  Its  radius 
is  Y centimeters.  The  susceptible  piece  part  is  located  a distance  of  Z centi- 
meters above  the  center  of  the  plane.  If  the  activity  of  the  dust  is  A photons/ 
cm2-sec.  and  only  a distance  attenuation  is  appropriate,  the  dose  rate  Dpp  at 
the  piece  part  is 


if 


» 1 


6 


pp 


For  a point  source  located  directly  beneath  the  piece  part,  the  dose  rate  is 


PP 


(73) 


The  correction  factor  to  convert  point  source  results  to  planar  source  results  is 


C.F. 


(7A) 


For  any  LRU,  then,  the  first  step  is  to  determine  the  critical  particle  size, 
and  the  pertinent  nondimensionalizing  constant.  From  the  figures  presented 
earlier,  a dose  at  1 meter  is  obtained.  This  dose  is  converted  to  the  dose 
(still  point  source)  at  the  piece  part  by  multiplying  by  the  factor  10‘*/Z^  . 
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This  dose  is  then  multiplied  by  the  correction  factor  shown  previously.  This 
dose,  then,  is  the  dose  which  is  accumulated  by  the  piece  part  of  a particular 
LRU  for  a given  TI,  aT,  etc.  This  dose  determination  is  done  for  all  the 
potentially  susceptible  piece  parts  of  the  LRU  in  question.  This  dose  for  each 
mission-critical  piece  part  is  added  to  the  inmersion  dose,  and  all  other  doses 
(including  any  possible  dose  resulting  from  exposure  to  prompt  radiation  during 
the  mission).  This  total  dose  is  compared  to  the  total  dose  susceptibility 
threshold  of  the  piece  part  (which  is  assumed  known).  If  the  dose  level  is 
below  the  susceptibility  threshold  for  each  mission-critical  piece  part,  then 
the  LRU  is  not  vulnerable  and  no  filter  is  required  to  protect  the  LRU  from 
total  dose  effects  from  the  accumulated  dose.  If  this  result  exceeds  the 
threshold  susceptibility  for  any  one  mission-critical  piece  part,  the  LRU  is 
vulnerable  and  protection  is  required  in  the  form  of  a filter. 

This  analysis  is  repeated  for  each  mission-critical  LRU.  The  most  vulnerable 
of  the  LRUs  will  drive  the  filter  selection.  At  this  point,  the  difference  or 
margin  between  the  susceptibility  threshold  level  and  the  accumulated  dose  level 
is  known  and  is  the  key  element  in  the  filter  criteria  selection.  This  margin 
is  related  to  the  amount  of  dust  which  must  be  trapped  by  a filter.  Since  in 
this  analysis  mass  and  dust  distributions  are  not  available,  the  approach  used 
in  the  cockpit  filter  criteria  selection  is  not  appropriate.  A different 
procedure  will  be  used. 

If  the  curves  depicting  the  mass  and  dose  distribution  with  particle  size 
for  the  cockpit  analysis  are  recalled,  it  is  apparent  that  the  LRU  and  cockpit 
act  as  perfect  filters  for  particles  whose  sizes  are  larger  than  the  critical 
particle  size,  R^.  Therefore,  the  perfect  filter  results  of  appendix  8 with 
the  appropriate  constant  (the  constant  in  this  case  must  be  based  on  the  mass 
rate  of  flow  of  air  to  the  LRU),  can  be  used  to  determine  the  necessary  filter 
trapping  efficiency.  The  procedure  is  to  convert  the  dose  margin  to  point 
source  (at-l-meter)  results  by  going  through  the  conversion  procedure  described 
above  in  reverse.  The  result  will  be  subtracted  from  the  dose  result  on  the 
pertinent  accumulated  dose  curve  at  the  10,000-micron  point  on  the  ordinate. 

This  result  corresponds  to  some  particle  size.  The  filter  must  trap  all  particles 
larger  than  this  size  to  provide  the  required  protection.  This  fixes  the  filter 
efficiency. 
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The  trapping  efficiency  of  the  filter  is  now  known.  The  dust  mass  and  the 
dust  mass  distribution  of  the  filter  now  can  be  found  from  the  perfect  filter 
results  of  appendix  B in  the  sane  manner  as  was  done  in  the  cockpit  analysis 
section.  Note  that  the  nondimensional  constant  in  this  case  must  be  based  on 
the  total  mass  rate  of  flow  of  the  cooling  air  through  the  filter  to  the  entire 
electronic  equipment  payload,  not  just  to  the  single  LRU. 

3.  EXAMPLE  LINE  REPLACEABLE  UNIT  ANALYSIS 

An  example  LRU  will  now  be  considered  as  an  illustration  to  clarify  the 
technique.  A case  of  T!  = 30  minutes,  aT  = 30  minutes  will  be  considered.  It 
is  assumed  that  the  LRU  of  interest  is  the  critical  one  in  the  aircraft,  based 
on  preliminary  analysis  of  the  type  described  previously.  This  example  was 
chosen  to  represent  a typical  LRU  of  one  of  today's  aircraft.  Its  physical 
configuration  is  presented  in  figure  22.  It  is  assumed  that  it  requires  a mass 
rate  of  flow  of  2 Ibm  (air)/minute.  Thus  from  equation  (A-44),  Kj  = 0.305, 
where  2 Ibm  (air)/minute  (converted  to  grams/hour)  replaces  m^^  in  the  Kj 
equation,  and  the  LRU  dimensions  replace  the  cockpit  dimensions.  The  constant 
has  been  multiplied  by  0.922  to  convert  it  to  rads(Si).  Applying  the  LRU 
parameters  to  equation  (31),  the  critical  particle  size,  R^  is  36p.  From 
figure  15,  the  dose  for  36u  at  1 meter  from  a point  source  is  30.5  rads(Si). 

For  this  LRU,  Z = 2 cm  and  Y = 8 cm.  Therefore,  the  dose  from  figure  15  is 
multiplied  by  the  factor  10‘*/Z2  to  correct  for  the  difference  in  at-l-meter  and 
actual  distance  calculations.  Next  multiply  by  the  correction  factor,  C.F.  = 
0.1733  (equation  (68))  to  account  for  the  fact  that  it  is  a distributed  source 
and  not  a point  source.  The  total  dose  then  which  impinges  on  the  susceptible 
component  is  1320  rads(Si).  It  is  assumed  that  the  piece  p rt  has  a suscepti- 
bility threshold  of  820  rads(Si).  Therefore,  the  margin  of  500  rads(Si)  must 
be  collected  by  the  filter. 

Now  to  determine  the  filter  requirements,  divide  the  margin  dose  of  500 
rads(Si)  by  the  two  correctioi  factors  above,  resulting  in  an  at-l-meter  point 
source  dose  of  1.15  rads(Si).  If  this  is  now  divided  by  Kj,  the  result  is  3.77, 
which  has  been  converted  back  to  rads(tissue)  and  then  nondimensionalized.  This 
value  can  now  be  used  directly  on  the  perfect  filter  graphs  in  appendix  B.  In 
this  case,  go  to  figures  B-35  and  B-44.  The  value,  on  the  figures  at  r = 10,000 
microns  for  AT  = 30  minutes,  is  24.  Subtracting,  24  - 3.77,  results  in  a level 
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of  20.23.  On  figure  B-35  for  aT  = 30  minutes,  this  level  occurs  at  a particle 
size  of  60  microns.  Therefore,  the  filter  efficiency  must  be  1 for  particles 
60  microns  or  larger,  i.e.,  the  filter  must  trap  all  particles  larger  than  60 
microns.  The  filter  dust  distribution  then  can  be  obtained  from  figure  B-7, 
using  a based  on  the  total  mass  rate  of  flow  required  for  all  electronics. 

Only  the  portion  of  the  distribution  above  60  microns  is  pertinent.  The  remain- 
der is  passed  by  the  filter.  The  total  dust  mass  in  the  filter  can  be  obtained 
from  figure  B-16,  the  dose  distribution  from  figure  B-26  and  the  total  dose 
accumulated  in  the  filter  from  figure  B-35  in  a similar  manner. 

4.  FILTER  LOCATION 

There  is  one  remaining  aspect  of  the  pr  jlem  which  must  be  considered  in 
this  section.  How  close  to  the  filter  (either  the  electronic  equipment  or  the 
cockpit  filter)  may  susceptible  electronic  equipment  be  located?  The  answer  to 
this  question  would  delineate  a “forbidden  region"  in  the  volume  adjacent  to  the 
filter  in  which  no  susceptible  electronics  may  be  located.  This  volume  would  be 
a strong  function  of  the  equipment  susceptibility  and  how  much  dose  it  accumulates 
from  other  sources. 

Consider  now  a filter  which  has  trapped  a considerable  amount  of  dust  during 
the  cloud  penetration.  The  filter  dose  at-l-meter  assuming  a point  source  of 
radiation  is  known  from  previous  considerations.  However,  for  nearby  electronic 
equipment,  the  dust  distribution  does  not  "look  like"  a point  source,  but  rather 
a distributed  source,  because  the  filter  element  is  assumed  to  be  a plane  per- 
pendicular to  the  flow  through  the  filter.  Therefore,  a correction  factor 
similar  to  the  one  developed  previously  for  the  LRU  plenum  chamber  is  needed. 

(It  is  noted  that  the  geometry  chosen  here  is  only  applicable  for  a planar-element 
filter.  Some  filters  have  other  geometries,  e.g.,  cylindrical  filter  elements, 
and  must  be  considered  separately,  but  in  a similar  manner.) 

This  geometry  is  indicated  in  figure  23  with  the  appropriate  distances 
shown.  The  piece  part  is  located  at  the  point  PP.  The  filter  element  is 
represented  by  the  circle  of  radius  d,  which  has  a uniform  dust  distribution 
and  specific  activity  A in  photons/cm^-sec.  The  differential  area  pdpdo  is 
located  at  point  DA.  The  dose  rate  at  the  piece  part  is 
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where 


= (a  + d)^  - 2p(a  + d)  cose  + 


(76) 


Performing  the  indicated  operations,  Dpp  is  obtained 


PP  " T a‘Ti'+  id) 


(a  + d 


■ 


(77) 


for  the  dose  rate  at  the  piece  part  due  to  the  distributed  dust.  If  all  the 
dose  were  concentrated  at  the  center  of  the  filter  element,  the  dose  rate  at 
the  piece  part  would  be 


6 

PP  ^ (a  t d)^ 

Therefore,  the  correction  factor  for  this  situation  is 

C.F.'  = (a 

d^  a(a  + 2d) 


(78) 


(79) 


The  dose  results  from  the  previous  work  must  first  be  corrected  for  the 
actual  distance  from  the  filter  element  to  the  piece  part  location,  or  any  loca- 
tion of  interest,  and  then  corrected  by  the  correction  factor  above.  Using  this 
procedure, contours  of  constant  dose,  or  isorads, may  be  determined.  The  equip- 
ment must  not  be  located  such  that  its  threshold  susceptibility  is  exceeded, 
i.e.,  inside  the  pertinent  isorad.  Generally,  this  forbidden  region  is  rela- 
tively small,  i.e.,  on  the  order  of  inches  away  from  the  filter,  because  of  the 
rapid  fall-off  of  dose  rate  with  distance  away  from  the  filter  and  the  relatively 
high  thresholds  of  the  equipment.  However,  this  forbidden  region  must  be  a 
consideration  in  equipment  and  filter  location  planning. 

The  determination  of  the  distance  which  the  filter  must  be  located  away  from 
the  crew  is  much  simpler.  Generally,  this  distance  is  at  least  a meter  or  greater 
and  at  this  distance  the  point  source  approximation  is  reasonably  accurate.  Only 
very  near  the  filter  must  the  correction  factor  above  be  used. 
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SECTION  VII 

MISCELLANEOUS  CONSIDERATIONS 

The  major  thrust  of  this  work  is  toward  the  protection  required  by  the  crevj 
and  the  electronic  equipment  from  the  contaminated  dust  ingested  into  the  air- 
craft via  the  environmental  control  system.  These  problems  were  addressed  in 
detail  in  the  preceding  sections.  This  section  will  briefly  address  the  aircraft 
engines  and  the  dose  due  to  dust  accumulated  on  the  aircraft  exterior. 

1.  ENGINE  INGESTION  OF  DUST 

The  air  flow  through  the  engines  is  much  larger  than  the  air  flow  to  the 
environmental  control  system.  Therefore,  much  greater  masses  of  dust  will  be 
ingested.  However,  because  the  engine  interior  is  streamlined  to  offer  low 
flow  resistance,  little  of  the  ingested  dust  will  be  accumulated.  There  is  no 
problem  of  the  engine  being  a source  of  radiation  due  to  trapped  dust.  The 
only  significant  potential  problem  is  that  the  large  amounts  of  dust  ingested 
could  erode  the  compressor-and-turbine  rotor  and  stator  blades,  and  cause 
unacceptable  damage. 

It  is  beyond  the  scope  of  this  work  to  investigate  any  potential  damage 
caused  by  the  dust.  All  that  will  be  done  is  to  present  the  technique  of 
determining  the  total  mass  and  the  mass  distribution  with  porticle  size  of  the 
dust  ingested  by  the  engines  during  the  cloud  penetrations.  It  is  hoped  that 
this  information  could  be  used  by  engine  design  and  test  engineers  to  ensure 
that  the  mission  completion  capability  of  the  engine  is  not  compromised  by 
cloud  penetrations. 

Since  the  engine  traps  little  or  no  dust,  all  of  the  dust  which  is  ingested 
flows  through  and  out  of  the  engine.  Recall  that  the  perfect  filter  analysis  in 
appendix  0 determined  all  of  the  dust  which  entered  the  filter  and  it  was 
assumed  to  trap  it  all.  Therefore,  the  perfect  filter  results  of  appendix  li 
are  pertinent  to  the  engine.  However,  in  this  case,  the  results  merely  show 
what  has  passed  through  the  engine  during  the  penetration  and  not  what  was 
trapped. 
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The  total  mass  ingested  by  the  engine  during  a penetration  can  be  obtained 
by  figure  B-23  of  appendix  B with  the  constant  (K^)  based  on  the  total  mass  rate 
of  flow  of  air  through  the  engine.  The  mass  distribution  of  the  dust  as  a func- 
tion of  particle  size  can  be  obtained  from  figures  B-5  through  B-22. 

2.  DUST  ACCUMULATED  ON  THE  AIRCRAFT  EXTERIOR 

The  dose  due  to  the  mass  of  dust  trapped  on  the  aircraft  exterior  is 
addressed  briefly  in  appendix  A.  Generally,  this  dose  is  relatively  small  and 
will  not  affect  either  the  crew  or  the  electronic  equipment.  The  primary 
interest  to  date  in  this  dust  is  related  to  the  recovery  of  the  aircraft  by 
maintenance  personnel  after  a mission  involving  a penetration  of  a rad’cactive 
dust  cloud.  Washing  the  aircraft  and  other  decontamination  techniques  should 
minimize  the  hazard  due  to  external  dust  accumulations.  However,  the  main 
problem  in  recovery  will  be  changing  the  filters  and  LRUs  with  open-cycle 
cooling.  These  problems  should  be  addressed  and  recovery  techniques 
formulated  to  minimize  the  hazard  to  ground  crews. 
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SECTION  VIII 
CONCLUSIONS 

The  purpose  of  this  report  was  to  investigate  the  penetratior  of  a raflio- 
active  cloud  by  a manned  aircraft.  The  concerns  associated  with  such  a penetra- 
tion are  three:  the  crew,  the  electronic  equipment,  and  the  engines.  If  any  of 

the  thre  • eyoerience  significant  performance  degradation,  the  strategic  mission 
of  the  aircraft  could  be  in  jeopardy.  Because  detailed  quantitative  results  are 
extremely  system  dependent,  specific  conclusions  about  the  effects  of  cloud 
penetration  on  crew  members,  electronics,  and  engines  cannot  be  drawn.  However, 
some  general  observations  can  be  made. 

1.  The  cloud  inmersion  dose  is  accumulated  equally  by  the  crew  and  the 
electronic  equipment  during  cloud  penetration.  The  only  feasible  action  to 
reduce  this  dose  is  to  avoid  the  cloud.  Shielding  to  attenuate  this  cloud 
inmersion  radiation  is  impractical  because  of  the  great  weight  of  material 
required  to  provide  any  appreciable  attenuation. 

2.  In  addition  to  the  cloud  immersion  dose,  the  crew  could  be  subjected  to 
potentially  severe  hazards  associated  with  the  dust  ingested  into  the  cockpit 

by  the  environmental  control  system  (ECS).  Although  these  hazards,  particularly 
the  ionizing  dose  accumulations  from  dust  accumulated  in  the  cockpit  and  skin 
burns  from  the  beta  radiation  associated  with  the  radioactive  dust,  can  be 
reduced  considerably  by  the  installation  of  a suitable  filter  in  the  ECS. 

3.  Although  precise  filter  criteria  for  the  cockpit  ECS  filter  are  the  end 
product  of  a detailed  analysis  ang  are  system  and  threat  dependent,  study  of  the 
cockpit  dose  results  in  appendix  D reve.ls  that  filters  which  trap  all  particles 
in  excess  of  6 to  8 microns  in  radius  would  provide  good  protection  to  the  crew. 
This  capability  is  well  within  the  present  state  of  the  art  in  filter  des  gn. 

4.  Although  the  ionizing  dose  accumulated  by  piece  parts  in  electronic 
equipment  is  system  and  threat  dependent,  it  is  observed  that  the  critical 
particl.e  size,  R^,  for  typical  plenum  chamber  geometries  and  mass  flow  rates 

of  cooling  air,  is  relatively  large,  i.e.,  50  microns  or  larger.  Because  fall- 
out rapidly  depletes  the  cloud  of  large  particles  and  the  plenum  chambers  of 
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representative  line  replaceable  units  trap  relatively  large  particles,  the 
amount  of  dust  accumulated  in  the  plenum  chambers  Is  relatively  small.  There- 
fore, the  piece-part  ionizing  dose  accumulation  is  probably  relatively  low. 
Therefore,  alleviation  of  this  potential  mission  crippler  should  be  relatively 
straightforward. 

This  report  has  shown  that  cloud  penetration  by  manned  aircraft  could  result 
in  unacceptable  performance  degradation.  Strategic  aircraft  should  be  subjected 
to  detailed  analyses  to  determine  potential  vulnerabilities  and  the  corrective 
actions  necessary  to  reduce  these  vulnerabilities. 
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APPENDIX  A 
UNIFORM  FALLOUT  MODEL 

This  appendix  develops  a preliminary  model  using  the  cloud  dust  density  and 
specific  activity  functions  which  were  presented  in  section  II.  The  relation- 
shM.s  from  section  II  are  independent  of  particle  size.  This  model  may  be  viewed 
as  a first  order  approximation.  Because  particle  size  is  not  a consideration, 
the  dust  cloud  effectively  may  be  modeled  by  a cloud  consisting  of  dust  particles 
which  are  uniform  in  size  and  activity.  This  model  does  take  fallout  into 
account,  but  in  a rather  simplified  fashion. 

For  the  purposes  of  this  development,  the  dust  density  and  specific  activity 
functions  are  written  in  the  form  developed  directly  by  curve  fitting  Whitaker's 
data  (section  II). 


Pd(t)  » a^(t*'"  + bf") 

(A-1) 

A(t)  = Ajf'‘ 

(A-2) 

where 


a.|  = 3.46  X 10"®  L^  gms(dust)/cm^ 
t = 1.2 
m = 1 .6 


n = 0.3 
b = 1 .313 


A = 4 22  X 10'^  photons 

■ gm(dust)-hour 
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1 . DUST  INGESTION  BY  THE  ENGINES 

The  mass  rate  of  flow  of  air  through  the  engines  Is  given  by  lii,.  grams  (air) 

o6 

per  hour  and  is  constant  for  a given  engine  and  flight  condition  of  the  aircraft. 
The  air  drawn  through  the  engines  contains  dust  in  the  ratio  Pj(t)/p  grams 

Q d 

of  dust  per  gram  of  air.  The  mass  rate  of  flow  of  dust  through  the  engines 
Mg{t)  in  grams  per  hour  is  given  by  the  relation 


f,  . i . ■«  *(  (,- . „-«) 

e'  ' ae  Pg  “^a  ' 

t^  < t £ t, 

(A-3) 

Me(t)  = 0 

all  other  t 

(A-4) 

For  some  purposes,  such  as  considerations  of  the  effect  of  engine  lifetime 
after  dust  ingestion,  one  msy  be  interested  in  the  total  mass  of  dust  ingested 
by  the  engines  during  a cloud  penetration.  This  can  be  obtained  by  Integration 
of  the  above  equation. 


Mg(t)  = 0 


•"ae  ®1 


Me(t)  = „ 


"’ae 

M,(t)  - 


t/-"’  - t 


srr 


,-m  /^,-n  . , )-n\ 

— ) 


fv-  - v-  _ (tr  - 
— spi ^ V — rs — /\ 


2.  DUST  ACCUMULATION  IN  AIR  CONDITIONING  FILTERS 


t < t.|  (A-5) 

t.j  1 t <.  t^  (A-6) 

t i.  t^  (A-7) 


Also  of  concern  to  this  analysis  is  the  mass  of  dust  collected  by  a perfect 
filter  in  the  bleed  air  line  from  the  engines  which  supplied  air  to  the  crew 
compartmcnt/electronics.  The  accumulation  of  radioactive  dust  in  the  filter  is 
a source  of  radiation  which  contributes  to  the  dose  received  by  the  aircrew 
members  and  electronic  equipment. 
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If  one  assumes  that  the  trapped  dust  Is  a point  source  of  radioactivity  and 
neglects  any  shielding,  one  can  calculate  the  dose  rate  and  dose  at  some  distance, 
dj  centimeters,  from  the  filter.  This  distance  represents  the  distance  to  the 
crewmember/electronic  equipment  from  the  filter. 

The  mass  of  dust  (Mf(t))  trapped  by  the  filter  can  be  calculated,  assuming 
100  percent  filter  trapping  efficiency,  in  a manner  similar  to  the  one  used 
above.  This  relation  is  given  by  the  following: 


Mf(t)  = 0 


m - 3i  j 

Mf(t)  = -5L_i 


rh  ^ da 
M.(t)  = -il-L 


* 1-m  *i-m  /*i-n  * i-n\ 


. i-m  . i-m  /.  i-ii  i-n\ 
‘i  ■ ‘f  . . Pf  - ‘i  1 
sn  ^ \ Tn  I 


t <. 


(A-8) 


< t < t^  (A-9) 


t >.  t^ 


(A-10) 


The  product  M^(t)  A(t)  is  the  effective  source  strength  of  the  filter  trapped 
dust.  The  dose  rate  measured  at  a distance  d centimeters  away  from  the  source  is 


Df(t)  = 


(A-11) 


where  C is  a conversion  factor  with  units  of  rads(tissu9)  centimeter^/photon. 
The  dost  rate  due  to  the  dust  trapped  in  the  filter  then  is  given  by 


i)f(t)  = 0 


tit. 


(A-12) 


i)f(t)  = 


Df(t)  = 


^ '"af  ^ 
4itd^ 


C A,  m^^  a. 


4„d2 


i-m.t'-"’ 

. i-m  , i-m  i-n  * i-n\ 


t^  <.  t s.  t^  (A-13) 


^ t >.  tf 


(A-14) 
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The  total  dof.'  accumulated  by  time  t due  to  the  filter  dust  is  obtained  by 
integration  of  the  dose  rate. 


Df{t)  = 0 t ;S.  t. 


C A,  m , a. 


where 
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3.  DOSE  DUE  TO  CLOUD  RADIOACTIVITY 

Another  primary  source  of  dose  to  the  crew  results  from  the  integrated 
effect  of  the  dose  from  the  radioactive  debris  surrounding  the  aircraft  as  it 
passes  through  the  cloud.  It  is  assumed  that  the  cloud  is  an  infinite,  home- 
geneous  mixture  with  a distributed  source  of  radiation  of  strength  given  by  the 
product  of  the  dust  density,  and  the  specific  activity  A(t).  For  simplicity, 
consider  the  aircraft  to  be  adequately  represented  by  a sphere  of  radius  R 
centimeters  and  neglect  any  shielding  effect  of  aircraft  structure.  For  any 
source  point  at  a distance  d > R from  the  center  of  the  sphere,  the  rate  at 
which  photons  are  emitted  from  a volume  element  dV  is  given  by  A{t)  dV. 

Since  the  source  is  assumed  isotropic  in  nature  A(t)  dV/And^  is  the  fluence 
at  the  center  of  the  sphere  from  the  volume  element.  One  must  also  account  for 
atmospheric  attenuation  of  these  photons  which  requires  a factor  of  exp 
■y-(u7o3)p3d^  where  u'/p^  is  the  mass  attenuation  coefficient  for  the  photons 
of  interest.  To  convert  from  photons  per  centimeter^  to  dose  units  use  the 
conversion  constant  C,  which  is  a function  of  the  energy  level  of  the  photons. 
The  total  dose  rate  due  to  all  volume  elements  is  given  by 

-Mpd 

f C P,(t)  A(t)e  4xd2  dd 

0 (t)  . J — ^ (A-1 

' ^ 4iid2 


or 


S(,(t) 


c p^(t)  A(t) 

u 


(A-20) 


since 


u ~ 10"5,  R ~ 10^  cm,  and  e”**  1 


therefore. 


C p.(t)  A(t) 

D (t)  =-^_, 

^ 11 


(A-21) 
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Subst'jtutinq,  the  equation  becomes 

6g(t)  =-^Aj  t'^  a,,  [f™  + bt'"]  t.  < t < tf 


D,(t)  = 0 


The  total  dose  accumulated  at  time,  t.  Is  given  by 

t 


all  other  t 


DJt)  = f 6 (t)  dt 
‘i 

Performing  the  indicated  integration 


D,(t) 

= 0 

t 

Oc(t) 

. *1  *i[m+FT  (*i 

P L 

l-(m+t)  _ ^l-(m+l)^ 

+ (t  i-(n+i) 

.i-(n+t)\l  t,  i t 

n+t-1  \‘i 

'•  /J  ' 

Oj.(t) 

^ ®i  r 1 f, 

y'  [_ra+l-l  V 

b i-(n+t)  i-Cn+trt] 
n+FT  Vi  ■ /J 


t >.  tf 


(A-22) 

(A-23) 


(A-24) 


(A-25) 


(A-26) 


4.  DOSE  DUE  TO  EXTERNAL  ACCUMULATION  OF  DUST 

If  radioactive  dust  is  accumulated  on  the  skin  of  the  aircraft  or  in 
crevices  and  discontinuities  caused  by  joints  in  the  external  structure,  it 
would  also  contribute  to  the  radiation  dose  received  by  the  crew/electronic 
equipment. 

Since  the  dose  from  this  source  would  be  Inversely  proportional  to  the  square 
of  the  distance  from  the  source  point  to  the  crew  nember/electronic  equipment,  it 
would  appear  reasonable  to  consider  only  the  dust  accumulated  in  the  near  vicin- 
ity of  the  crew  member/electronic  equipment.  For  example,  in  determining  the 
dose  to  the  pilot,  one  would  restrict  bis  attention  to  the  accumulation  of  dust 
in  the  cockpit/crew  station  area. 
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If  one  could  reasonably  estimate  the  amount  of  dust  accumulated  in  the 
pilot's  vicinity  as  G grams  at  an  effective  distance  of  d centimeters  one  could 
estimate  the  dose  due  to  the  source.*  To  be  conservative  one  might  assume  that 
all  of  this  dust  was  accumulated  at  the  time  of  entry,  t^.,  into  the  radiation 
cloud  and  remained  there  during  the  remainder  of  the  mission. 

Then  the  dose  rate  at  position  R from  this  source  would  be  given  by 


D (t)  = 

4nd^ 


-i 

Da-(t)  = ^t 

4nd^ 


The  integrated  dose  then  is 


t i.  t. 


Dae(t)  = 0 


GCA, 


t <.  t. 


t >.  tj 


(A-27) 


(A-28) 


(A-29) 

(A-30) 


This  area  has  been  investigated  experimentally  and  analytically  by  several 
researchers.  No  attempt  will  be  made  here  to  correlate  the  experimental  find- 
ings with  this  rather  crude  first  cut  analysis.  Future  effort  in  this  area  is 
needed.  However,  in  this  work  this  dose  is  much  less  critical  than  the  filter 
and  cloud  doses,  and  has  been  presented  briefly  only  for  the  purposes  of 
completeness. 

5.  TOTAL  DOSE 

The  total  radiation  dose  due  to  penetration  of  the  cloud  then  is  the  sum  of 
individual  doses. 


*As  a "first  cut"  at  the  accumulation  of  dust  on  the  aircraft  exterior,  one 
could  measure  or  estimate  the  volume  of  all  the  cracks,  crevices,  etc.,  and  the 
volume  of  surface  accumulation,  which  would  be  a function  of  the  boundary  layer 
thickness.  If  then  the  dust  density  in  the  cracks,  etc.,  is  known  or  assumed, 
a total  mass  can  be  approximated. 
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DT-(t)  = Df(t)  + D^.(t)  + Dgg(t) 


(A-31) 


6.  SPECIFIC  RESULTS 

Substituting  the  values  for  the  cloud  parameters,  i.e.,  m,  n,  C,  and  b,  into 
the  equations  developed  previously,  and  grouping  the  remaining  constants,  the 
following  equations  are  obtained. 


^t) 

S/2 


• = + 1.313  t-°*3 


< t < t^  (A-32) 


(t.-“-®  - t-“*s  1 ,,,  , 

e/T  - V 1 / 1 1.313  Un.i  t ft.7\ 

tj  <.  t £ t^ 

(A-33) 

^IfZ 

le/f(t)  _ (t^-®*®  • V®-®)  ^ 1 313  , . 

-u  ^ j 

t >.  t.f. 

(A-34) 

'“1/2 

D t->-2 

(A-35) 

t^  < t <.  t 

^ ■ [sT  (V"-‘ - V'V  ^ 

■’)].-■*  tit, 

(A-36) 

Df(t) 

S 


6.25  t.""'®  + 2.083  f®-®  + 3.751  t®-®  + 9.379  t^®-^  f®*^ 

8.33  t^-°-®  t-®-2  - 13.13  t^®-®  < t < f,. 


(A-37) 
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D.(t) 


= 6.25  + 13.13  - 6.25  - 13.13 


+ (9.379  + 8.333  - 8.333  - 9.379 


D.(t)  , \ 

c . = (t-*-®  + 1.313  f*- 


K4 


DJt) 


1 /*.  -1.8  ^-i.e\  , 1 .313  -0.5  X 


Djt) 


' , 1 /t  -1.8  ..  -;.8\  , 1.313  -0.5  + -o.5\ 

“ ■ 13  ) * “O-  1*1  - *f  ) 


K, 


Ml_l= 


D=o(t) 


■ = 5 


where 


K = '”ae  ^1  gm(dust) 

* p nr 

d 

y _ '"af  gm(dust) 

''■2  ■ p hr 

d 

_ C Aj  a^  rads(tlssue) 
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^ rads(tissue) 

■ p' 

= ^-XaMtlssuel 
4Tid^ 

Note  that  the  expression  for  the  amount  of  dust  ingested  into  the  engine  is 
functionally  identical  to  the  expression  for  the  amount  of  dust  trapped  by  the 
perfect  filter.  The  only  difierence  is  in  the  constant,  i.e.,  Kj  corresponds  to 
the  dust  ingested  by  the  engines  and  to  the  dust  trapped  by  the  filter. 
Therefore,  these  two  equations  were  combined  in  equations  A-32  through  A-34. 

7.  GRAPHICAL  RESULTS 

These  equations  are  all  completely  analytical  and  easy  to  solve  for  any 
given  set  of  input  parameters.  Representative  results  have  been  obtained  and 
are  shown  in  figures  A1  through  AlO.  The  ordinate  in  each  case  is  the  inde- 
pendent variable  of  interest  nondimensionalized  by  the  pertinent  constant,  i.e., 
Kj,  Kj,  Kj,  K^,  or  Kj.  These  results  then  are  general.  For  a specific  aircraft, 
threat,  and  mission  profile,  the  constants  can  be  determined  and  the  specific 
results  can  be  determined. 

Figure  A1  depicts  the  mass  rate  of  flow  of  the  dust  to  the  engines  or  the 
filter.  Fallout  of  the  dust  in  the  cloud  is  reflected  in  this  figure  through 
the  decrease  with  time  of  the  mass  rate  of  flow  of  the  dust.  For  a no-fallout 
situation,  the  mass  rate  of  flow  of  dust  would  be  constant.  Therefore,  fallout 
is  a significant  factor  in  this  development.  Figure  A2  depicts  the  total  mass 
ingested  by  the  engines/trapped  by  the  filter  for  two  cloud  entry  times  (TI), 

10  minutes  and  60  minutes  after  detonation.  The  total  mass  ingested/ trapped  as 
a result  of  a particular  penetration  duration  (aT)  is  read  from  the  ordinate  at 
an  exit  time  of  TF  = TI  + aT,  for  the  curve  corresponding  \o  the  proper  TI. 

Figures  A3  and  A4  depict  the  dose  rates  (at  1 meter)  due  to  the  dust 
trapped  in  the  filter.  Figure  A3  corresponds  to  a TI  of  10  minutes  and  figure 
A4  a TI  of  60  minutes.  Each  figure  depicts  a family  of  curves  corresponding  to 
penetration  durations  of  2,  5,  10,  15,  30,  45,  60,  120,  300,  and  600  minutes. 

The  rapid  decrease  in  the  dose  rates  with  time  is  due  to  two  factors,  the  fellout 
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of  the  dust  particles  from  the  cloud,  and  the  decrease  in  photon  emission  due  to 
the  normal  decay  process.  Note  that  the  dose  rate  does  not  go  to  zero  after 
exit  from  the  cloud  because  the  dust  trapped  in  the  filter  remains  in  the  filter 
after  cloud  exit  and  is  a source  of  ionizing  dose. 

Figures  A5  and  A6  depict  the  ionizing  dose  (at  1 meter)  due  to  the  dust 
trapped  in  the  filter  for  cloud  entry  times  of  10  and  60  minutes  after  detona- 
tion, respectively.  The  curves  in  the  figure  correspond  tc  the  same  penetration 
durations  given  previously.  Note  that  the  dose  rate  and  "^ose  results  are 
presented  for  a distance  from  the  filter  of  1 meter.  These  results  can  be 
related  to  other  distances  by  multiplying  the  results  by  tlie  factor  (1  meter)^/ 
(distance  in  meters)^.  Figure  A7  depicts  the  dose  rate  due  to  the  aircraft 
being  immersed  in  the  cloud.  (The  crew  and  all  aircraft  components  are  exposed 
equally  to  the  photons  emitted  from  the  cloud,  because  the  shielding  provided 
by  the  aircraft  is  negligible.) 

Figure  A8  depicts  the  dose  accumulated  by  the  crew  due  to  cloud  inmersion 
for  TIs  of  10  and  60  minutes  after  detonation.  The  dose  accumulation  is  read 
by  choosing  the  curve  corresponding  to  the  TI  of  interest,  determining  the  cloud 
exit  time,  TF  = TI  + aT,  and  reading  the  dose  from  the  ordinate.  Figure  A9 
depicts  the  dose  rate  (at  1 meter)  per  gram  of  dust  accumulated  on  the  aircraft 
exterior.  Figure  A10  depicts  the  dose  (at  1 meter)  per  gram  of  dust  accumulated 
on  the  aircraft  exterior  for  TIs  of  10  and  60  minutes  after  detonation. 

This  model  is  simple,  easy  to  use,  and  its  associated  equations  are  completely 
analytical.  For  any  given  set  of  input  parameters,  i.e.,  Ks,  TI,  and  oT,  results 
can  be  obtained  with  only  a slide  rule.  However,  this  model  does  not  realistic- 
ally represent  the  actual  radioactive  dust  environment  and  its  behavior  with  the 
time  and  may  not  yield  realistic  results. 
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Figure  A3.  Dose  Rate  from  Dust  in  Filter,  TI  = 10  Minutes 
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APPENDIX  B 

IMPROVED  FALLOUT  MODEL 

In  the  previous  appendix,  a first  cut,  simple  dust  cloud  model  was  Jeveloped. 
Because  this  model  does  not  account  for  dependence  of  the  results  on  particle 
size,  the  predicted  results  will  not  support  additional  detailed  analyses.  The 
lack  of  particle  size  dependence,  in  particular,  handicaps  the  filter  (.rapping 
analysis.  For  these  ic'^ons,  it  was  decided  to  develop  a more  realistic  model. 
This  model  and  its  predictec  results  will  be  investigated  to  determine  if  its 
more  realistic  results  (compared  to  the  Uniform  Fallout  Model)  offset  the  added 
complexity  of  the  model  and  associated  equations.  Although  the  equations  asso- 
ciated with  this  approach  appear  to  be  rather  complicated,  in  reality,  they  are 
relatively  simple  and  straightforward.  Relatively  simple  numerical  techniques 
are  employed  to  solve  them  using  only  a minute  or  so  of  computer  time,  rather 
than  the  hours  required  for  the  complex  dust  environment  codes. 

It  is  well  known  that  very  small  particles,  i.e.,  r < ly,  tend  tt  lain 
suspended  in  the  atmosphere  almost  indefinitely  (ref.  9).  One  method  ..i-  ultimate 
removal  is  by  their  being  used  as  precipitation  centers  for  water  vapor  and  fall- 
ing out  with  rain  or  other  precipitation.  The  point  is  that  they  do  net  fall 
out,  at  least  over  the  time  period  of  interest.  The  settling  veloci^i  of  very 
small  particles  is  expressed  in  Stokes  (ref.  10)  as 

= 1.2  X 105  r^  (E.-,) 

where  is  the  settling  velocity  in  cm/sec,  is  the  density  of  the  particle 
in  gm/cm^  and  r^  the  particle  radius  in  centimeters.  Substituting  in  the  value 
of  for  a representative  soil,  Pp  “ 2 gm/cm^  and  converting  the  radius  to 
microns,  then 


Vj  = 0.024  (B-2) 

where  r is  the  particle  radius  in  microns.  The  equation  was  derived  with  the 
assumption  that  the  Reynold's  Number  based  on  the  particle  diameter  and  settling 
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velocity  was  very  small,  i.e., 


2 X radius  x V 

Re  = =■  « 1 (B-3) 

''air 

where  is  the  viscosity  of  air.  The  relation  is  very  accurate  for  the  range 
of  particle  radius  from  1 micron  to  50  microns  (ref.  9).  Although  Stake's  Law 
was  derived  for  Re  « 1 , it  has  been  found  experimentally  that  there  is  good 
agreement  between  the  theory  and  experiment  up  to  RE  ~ 1 , where  divergence 
begins,  but  agreement  is  still  adequate  up  to  RE  ~10  (ref.  11).  Therefore,  the 
settling  velocity  concept  is  reasonably  accurate  up  to  a particle  radius  of  about 
lO^y,*  Based  on  review  of  available  literature,  the  initial  particle  size 
distribution  was  assumed  to  be 


f ~r"^*^  (B-4) 

with  cutoffs  at  r = O.lu  at  the  lower  end  and  r = 10,000u  at  the  upper  end. 

These  cutoffs  were  based  on  intuition  sharpened  by  study  of  available  information 
The  upper  limit  was  suggested  by  Whitaker's  work  (ref.  1)  and  the  lower  based  on 
consideration  of  possible  condensation  patterns  after  a detonation.  The  assump- 
tion is  made  that 


f(r)  = N r‘^*^  (B-5) 

where  f(r)  is  the  total  number  of  particles  with  radius  r in  a cubic  centimeter, 
and  N is  the  constant  of  porportional ity.  Then  the  mass  per  cubic  centimeter  of 
all  the  particles  with  radius  r is 


I"’ 

cm 


/V  4 fc,  -0.5 
Pd(r)  = N-rr 


(B-6) 


*Although  this  concept  breaks  down  for  particle  sizes  greater  than  lO^y,  equation 
(B-2)  was  used  in  the  numerical  integrations  used  to  calculate  the  results  in 
this  appendix.  Therefore,  the  results  associated  with  particle  sizes  in  excess 
of  lOOOii  are  somewhat  inaccurate. 
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where  is  the  density  of  each  solid  particle  (assumed  to  be  ?.  gm/cm^)  and 
is  the  standard  dust  density  distribution  (equation  (3)  in  section  II). 

It  has  been  shown  that  the  larger  particles  tend  to  fall  out  faster  than  the 
small  particles  (ref.  9).  It  is  assumed  that  all  particles  with  the  largest 
radius  fall  out  before  any  of  the  next  smaller  radius.  This  is  referred  to  as 
the  "nibbling  mouse"  assumption,  because  the  only  change  to  the  number  distribu- 
tion or  the  density  with  time  is  the  disappearing  end  point.  The  end  point  can 
be  envisioned  to  be  in  the.  process  of  hejng  eaten  by  a mouse  nibbling  the  curve 
at  a "rate,"  R(t).  Figures  B1  and  B2  show  the  number  distribution  and  the 
density  with  cutoffs  noted.  These  curves  were  obtained  by  fixing  the  constants 
involved  in  the  following  way.  The  following  relation  is  forced  to  hold. 

R(t) 

pj(t)  = P(j(r)  dr  (B-7) 

where  p^  is  defined  by  equation  (3)  in  section  II.  Note  that  the  time  dependence 
of  the  density  distribution  with  r enters  only  through  the  upper  limit  of  the 
integral.  The  value  of  R(t)  at  10  minutes  (or  0.167  hour)  is  defined  to  be  10'* 
microns.  The  above  relation  evaluated  at  10  minutes  after  detonation  fixes  the 
value  of  the  constant  N.  Therefore, 

f(r,t)  = 1.2  X 10“>  a^.  r'**®  0.1  <.  r <.  R(t)  (B-8) 

Pj(r,t)  = 9.9  X 10'^  a,  r'"*®  = 9.9  X 10'^  0.1  < r <.  R(t)  (B-9) 

where  K.^  = a.j  gm(dust)/cm3,  a nondimensional izing  constant.  Recall  that 

R(t) 

pj(t)  = 9.9  X 10'^  a^  y r'°-^  dr 
0.1 
or 

Pj(t)  = 0.2  a.  [R(t)°-^  - 0.3ie]  (B-10) 

To  determine  R(t)  note  that  equation  (B-10)  must  hold  for  all  t.  Therefore, 
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a.  + 1.313  r“'^)  = 0.2  a^  [R(t)’/^  - 0.316]  (B-11) 

which  results  in 

R(t)  = [s.04  + 1.313  + 0.316]^  (B-12) 

The  density  distribution  function  is  completely  defined  as  a function  of 
particle  size  and  time;  the  specific  activity  must  now  be  determined.  It  appears 
reasonable  to  assume  that  the  activity  of  the  smaller  particles  is  proportional 
to  the  entire  particle  mass  which  is  predominately  condensed  bomb  debris,  and 
that  of  larger  particles  is  proportional  to  the  mass  of  the  surface  coating  on 
the  particle  which  is  condensed  bomb  debris.  If  this  assumption  is  made  then  the 
specific  activity  per  particle  is  given  by 

Vticle-'-'  0.1<r<5  (B-13) 

Vrticle'-'-^  6<r<10^  (B-14) 

where  5 is  the  somewhat  arbitrary  boundary  (assumed  from  hereon  to  be  20p) 
separating  the  region  where  particles  are  assumed  to  be  composed  entirely  of 
condensed  debris  and  the  region  where  particles  are  assumed  to  be  coated  with 
bomb  debris. 

Relating  the  activity  to  the  mass  of  each  particle,  then 

A(r,t)  = fj  t‘‘-^  0.1  < r < 6 (E-15) 

A(r,t)  = fj  r*'  t‘‘'^  airslO"  (B-16) 

Making  use  of  the  condition  that  at  r = «, 

A(r,t)  _ = A(r,t)  ^ (B-17) 

r s r ->•  S 

This  results  in  the  relation 

fj  = fj«  (B-18) 
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To  evaluate  the  constant,  one  other  condition  Is  needed.  The  condition 
chosen  was  based  on  the  physical  situation.  It  is  assumed  that  the  10  minute 
conditions  in  all  the  fallout  models  considered  are  Identical.  Therefore,  the 
total  number  of  fission  fragments,  the  photons  per  unit  area  of  an  iiimerseri 
body,  and  the  dose  rates  for  both  cases  are  equal  ^ 10  minutes  after  detonation. 
For  the  Uniform  Fallout  Model, 

Dc(tl  = 4pd<*^  (B-19) 

For  the  Improved  Fallout  Model 


i 


D (t)  = 

C I 

u 


y 0(j(r)  A(r,t)  dr 
LO  1 


R(t) 

y Pj(r)  A(r,t)  dr 
6 


and  it  was  assumed  that 


(B-20) 


► 


0 (0.167)  = D (0.167)  (B-21) 

Improved  uniform 

Performing  the  necessary  operations  and  simplifying,  the  following  is  obtained. 


99.684, 

f,  i 11.83  A, 

‘ 2d'>-  - 0.01  « - 0.316 


f^  = 6 a = 236.55  A^ 


where  Aj  is  defined  in  appendix  A and  4 = 20y.  Therefore, 

A(r,t)  = 11.83  Aj  t'*'^  0.1<r<20u 

236.55  A,  t'*'^ 

A(r,t)  20u  <.  r <.  10'*)i 


(B-22) 

(3-23) 


(B-24) 

(B-25) 


Note  that  the  specific  activity  is  applicable  from  r = 0.1  to  r = lO**  for  all 
t.  The  "nibbling  mouse"  is  not  applied  to  the  activity  because  the  activity  is 
per  gm  of  dust. 


116 


AFWL-TR-73-82 


The  cloud  parameters  are  now  defined  and  using  the  same  assumptions  mentioned 
earlier,  i.e.,  zero  wind,  etc.,  the  mass  rate,  mass  dose  rate,  and  dose  calcula- 
tions performed  earlier  are  repeated  using  the  newly  derived  cloud  characteris- 
tics, which  are  now  particle  size  dependent. 

The  mass  trapped  in  the  filter  as  a function  of  tine  and  particle  size  is 
determined  assuming  that  all  the  dust  is  trapped  by  the  filter. 

Mf(r,t)  = i^^Pj(r,t)  (B-26) 


-2  -0  5 

= 9.946  X 10  ^ r 


O.lg  i r <.  R(t) 


(B-27) 


where  is  the  same  as  defined  in  equation  (A-44)  in  appendix  A. 

The  total  mass  trapped  is  the  integral  over  time  of  the  above  relation. 

Mx(r,t)  . . 

= 9.946  X 10"  r‘  -Mt  - ti)  0.1yir<R(t)  (B-28) 


Mf(r,t) 


■g—  = 9.946  X 10'^  (t„  - t,)  R(t)  < r^  < RCt,-)  (B-29) 


where  t^  <.  t^  <.  t,  and  r^  = R(tn).  Equation  (B-29)  is  handled  with  numerical 
techniques  for  plotting  purposes. 

This  combination  of  analytical  and  numerical  techniques  is  necessary  because 
of  the  "nibbling  mouse."  Figure  B3  is  representative  of  the  mass  rate  of  flow 
of  dust  into  the  filter,  M^{r,t).  Note  that  it  is  a simple  prism  with  the 
"nibbling  mouse"  eating  away  at  the  edge.  Figure  B4  is  a representation  of  the 
mass  of  dust  collected  in  the  filter.  This  three-dimensional  figure  is  consid- 
erably more  complex  than  the  previous  figure.  "Slices"  of  three-dimensional 
figure  B4  for  various  (t;.  tf)  sets  are  presented  in  figures  B5  through  B13. 

This  M^(r,t)  equation  set  can  be  integrated  over  r and  then  evaluated  at  t^ 
to  give  an  equation  set,  which  represents  the  cumulative  mass  (CMf)  as  a 
function  of  r,  for  some  (t,,  tf). 
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—A = 0.1989  (r“«  = 

'■2  ' 


0.3162) (tf-t^) 


0.1<r^R(tf)  (B-30) 


CM^(r)  r c 1 

= 0.1989  [R(tf)°-^  - 0.3162J  (tf-t^) 

R^i ) 

+ 9.946  x 10“^  X r r R(tf ) < r <.  R(ti ) (B-31 ) 

RUV)  " " 

Note  that  the  second  equation  must  be  integrated  numerically.  These  results  are 
presented  in  figures  B14  through  B22.  Figure  B23  presents  the  total  mass 
collected  by  the  filter  as  a function  of  time  after  weapon  detonation  for  the 
TIs  of  interest.  The  total  mass  results  In  this  figure  are  identical  to  the 
results  in  figure  A2  in  appendix  A. 

These  results  may  be  used  to  directly  determine  the  total  mass  (and  its  size 
distribution)  ingested  by  the  engine  in  the  same  manner  as  discussed  in  the 
Uniform  Fallout  I'odel  Analysis. 

Now  considei  the  dose  rate  and  dose  which  are  associated  with  the  dust 
trapped  in  the  filter.  This  dust  acts  as  a source  of  photons  which  could  affect 
the  crew/alectronic  equipment.  The  dose  rate  can  be  written  as 

Df(r.t)  = -l-H.(r.t)  A(r,t)  (B-32) 

Substituting,  the  following  set  of  equations  is  obtained. 


t^  < t i t^ 


Df(r.t) 


Df(r.t) 


Df(r,t) 


1.176  r'°-^  (t  - to  t'*'^' 
23.52  r'*-®  (t  - ti)  t''-^ 
23.52  (tn  - t^) 


0.1  < r <.  20 
20  <.  r < R(t) 
R(t)  < r^  < R(ti) 


(B-33) 

(B-34) 


{B-3E) 
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t i t- 


0 5 15 

= 1.176  r (tf  - t-i)  t'^'^ 


Df(r.t) 


23.52  r"^*^  (tf  - ti)  t‘ 


0.1  < r < 20 


(B-36) 


20  < r < R(tf)  (B-37) 


Df(l^jt)  -15  12 

-4 23.52  (tn  - t^}  t 


R(tf)  < < R(ti)  (B-38) 


From  these,  the  dose  distribution  as  a function  of  size  r for  a (t,-,  tf)  set 
can  be  obtained  by  integrating  over  t.  The  distribution  can  be  determined  at 
any  time  after  detonation.  The  time  of  most  interest  here  is  the  time  of  mission 
completion  (T^^)-  complete  set  of  equations  for  t^^  > t^  and  R(tf)  > 20u  is 

0.1  <,  r < 20 


D^{r) 


20  i.  r i R{tf) 


D.(r) 


V=  117.6  r->-=  [l.25  (t/-«  - t.”-®)  + (ti  - tf)]  (B-40) 


R(tf)  <.  r^^  i R(ti) 


-1^=  117.6  r;'-^  [l.25  (t;-«  - (t,  - t,)]  (B-41) 

where  r^  = R(tn)  and  t.  <.  t^  i t^. 


119 


1 


AFWL-TR-73-82 


If  R(tf)  < 20y,  then  the  above  set  of  equations  becomes 
0.1  < r < R(tf) 


5.88  r‘“-^  [l.25  (t^"*®  - t^“*®)  + (ti  - tf)  (B-42) 


R(tf)  :!  r < 20 


^ = 5.88  [l  .25  (t/-®  - t.''’®)  + (t,  - t„)]  (B-43) 


20  <.  r < R(ti) 


117.6  [l.25  (t;*^  - (t,  - tn)]  (B-44) 

These  equations  are  plotted  in  figures  B24  through  B32. 

Integrating  this  set  over  r gives  a cumulative-dose-as-a-function-of-r  set 
of  equations.  Note  that  again  numerical  integration  must  be  used  in  parts  of 
the  evaluation.  For  R(tf)  >.  20u, 

0.1  < r <.  20 


CD,  (r) 


= 11.76  [l. 25  (t/-®  - t.“-®)  + (ti-tf)]  (r“-®-0.316)  (B-45) 


20  i.  r <.  R{tf) 


CD,  (r)  CDf  (20) 


+ 235  1^1.25  (tf"*®  - t^"*®) 

^ \k°'^  (ti-tf)]  [o.2236  - r'"*®] 


120 


AFWL-TR-73-82 


R(tf)  <■  <.  R(ti) 


where  = R{tp)  and  <.  <. 

If  R(tf)  s.  20u,  the  preceding  set  of  equations  becomes 

0.1  i r <.  R{tf) 


(B-47) 


CDf  (r) 


11 


.76  [l.25  (t^°‘®  - t/-®)  + (ti-tf)]  (r°‘®-0.316)  (B-48) 


R(tf)  i r < 20 


CO, 


^(r)  CDf^^R(tf)) 


20 


Kj  K3 


+ 5.88  S fl.25  (t  ®-®  - 

R(tf)  L ^ ’ 


+ 1 


me 


(ti-tn)]  or 


(B-49) 


20u  i r <.  R{t.j ) 


CD,  (r)  CD,  (20)  R(ti) 


+ V'’'' 


(8-50) 


where  r^  = R(tn)  ~ equations  are  plotted  in  figures  B33 

through  B41  in  semilog  scaling  for  better  accuracy  at  the  higher  dose  levels  and 
figures  B42  through  B50  in  log-log  scaling  for  better  accuracy  for  the  lower 
dose  levels. 
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Another  useful  set  of  equations  can  be  determined  by  taking  equations  (B-33) 

' through  (B-38)  and  integrating  out  the  size  dependence  to  obtain  the  dose  as  a 

function  of  time.  Again  a numerical  technique  must  be  used,  and  r^^  = R(tp) 
with  t^.  i £ t^.  For  t.  <.  t i t^  and  if  R(t)  > 20p 

Df  (t) 

= (l01  .46  - 235.27  R(t)'“‘®)(o.25  t®‘®  - 1.25  t^°'®  + t.f°-^) 

I + 117. r (t^'®'^  - t'“-^)£  r„'*-®  (tn-ti)  Ar  (B-51) 

I R{t) 

if  R(t)  < 20y 


) 


D,  (t) 


= (n.76  R(t)“‘^  - 3.718s)(o.25  t“’®  - 1.25  t^“‘®  + 


20 


R(t) 

R(ti) 

+ 117.6  E (tn-ti) 

20 


ar 


(B-52) 


For  times  after  t^,  i.e.,  t > t^,  the  equations  become,  for  R(tf)  > 20y, 


^ R(ti) 


(tf  - t^)  (101.46  - 235.27  R(tf)'°'^) 


117.6  Z ^ r (tn  - t^-) 
R(tf) 


ar 


Df  (tf) 


(B-53) 


If  R(tf)  < 20u,  this  set  of  equations  becomes 


0,(t) 


20 


(tf  - ti)(ll.76  R(tf)“‘®  - 3.7185) 


■:  5.88  Z r.  (tn  - t^)  ar 
R(tf)  " 

R(ti) 

+ 117.6  Z fn'‘''  - tl) 


ar 


20 


Df  (tf) 


(B-54) 


L 
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The  filter  dose  as  a function  of  time  is  plotted  for  several  TI  and  TF 
combinations  in  figures  B51  through  B59. 

The  dose  rate  in  a cloud  is 

A(r,t)  (B-55) 

or 

-0  5-12  , V 

— = 1.176  r t 0.1o<r<20  (B-56) 

D.(r,t) 

-2-p = 23.527  r t 20  < r < R(t)  (B-57) 

\ 

Integrating  over  r 

47.054  t'**^  (o.4314  - R(t)"'''^)  (B-58) 


The  dose  then  is 


47.054 


2.157  it 


(B-59) 


The  dose  rate  and  the  dose  are  presented  in  figures  B60  and  B61 . 

In  all  these  relations  the  sunmation  symbol  indicates  numerical  integration, 
and  the  integral  symbol  analytical  integration.  Extensive  checks  of  the 
numerical  techniques  used  were  made  to  ensure  ••easonable  accuracy.  Comparing 
figures  B33  and  B51 , it  is  seen  that  the  two  different  approaches  yield  results 
at  t = 30  hours  very  close  in  value,  i.e.,  to  about  1 percent. 

It  is  noted  again  that  these  results  are  applicable  for  the  100-percent 
efficient  filter. 
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GRAPHICAL  RESULTS 

The  equations  associated  with  this  model  are  considerably  more  complex  than 
those  associated  with  the  simpler  Uniform  Fallout  Model,  and  dictate  the  use  of 
numerical  techniques  to  obtain  numerical  results.  Therefore,  extensive  graphical 
results  are  presented  to  support  accurate  interpolation.  This  model  is  more 
realistic  than  the  Uniform  Fallout  Model,  but  could  be  further  refined  by  inclu- 
sion of  the  dependence  of  the  dust  density  on  altitude  as  a function  of  time  and 
the  use  of  a better,  more  mathematically  precise  number  distribution,  i.e., 
refine  the  "nibbling  mouse"  assumption.  The  above  refinements  might  provide 
more  accurate  results,  but  the  increase  in  accuracy  would  probably  not  be 
significant. 

The  graphical  results  themselves  will  now  be  discussed  to  illustrate  their 
use.  Recall  that  the  filter  used  in  this  development  was  a perfect  filter.  It 
collects  all  ingested  dust  particles  for  the  penetration  Al . When  the  aircraft 
exits  the  cloud  at  TF  = TI  + iT,  no  more  particles  are  collected.  The  trapped 
mass  distribution  results  are  shown  in  figures  B5  through  B13.  Each  graph 
presents  ;he  dust  mass  collected  by  the  filter  as  a function  of  particle  size 
for  a particular  TI  and  with  penetration  duration,  iT,  as  a parameter.  Each 
graph  considers  aTs  of  2,  5,  10,  15,  30,  60,  120,  and  300  minutes.  These  fig- 
ures are  for  cloud  entry  times  ranging  fron  10  minutes  to  5 hours  after  weapon 
detonation,  which  should  cover  the  range  of  interest  for  most  analyses.  Particle 
size  distributions  range  from  0.1  micron  to  10,000  microns  radius.  A comparison 
of  figure  B5  for  a cloud  entry  time  of  10  minutes  with  figure  B13  for  a cloud 
entry  time  of  5 hours  shows  the  appreciable  effect  of  fallout  on  the  results. 

At  TI  = 10  minutes  a sizable  portion  of  the  .....ss  consists  of  particles  with 
r > 10  microns;  whereas,  for  TI  = 5 hours  all  of  these  particles  have  fallen  out 
prior  to  aircraft  entry,  and  thus  no  particles  in  this  size  range  are  collected 
by  the  filter. 

The  second  series  of  perfect  filter  graphs,  figures  B14  through  B22,  presents 
the  cumulative  mass  collected  by  the  filter  for  a particular  TI,  as  a func  in 
of  particle  size  and  in  terms  of  the  same  range  of  parameters  AT.  CurnuU-tive 
mass  CM^(r),  means  the  dust  mass  collected  by  the  filter  represented  by  dust 
particles  which  have  sizes  less  than  or  equal  to  r.  If  M^(r)  represents  the  dust 
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mass  distribution  function  for  the  filter  as  a function  of  particle  size  r for 
a TI,  aT  case,  then 

r 

CMf(r)  = f M.(r)  dr  (B-60) 

Figures  B14  through  B22  present  these  results  for  cloud  entry  times  ranging 
from  10  minutes  to  5 hours.  To  further  clarify  the  meaning  of  these  figures, 
consider  for  a moment  figure  B14.  This  figure  gives  the  cumulative  mass  collected 
by  the  filter  as  a function  of  particle  size  for  a cloud  entry  time  of  TI  = 10 
minutes.  For  illustration  purposes  select  a penetration  duration  of  aT  = 30 
minutes,  then  focus  on  the  question  of  how  much  dust  will  be  collected  which  has 
particle  sizes  ranging  up  to  50  microns.  Entering  the  abscissa  at  50  microns  and 
moving  up  to  the  curve  corresponding  to  AT  = 30  minutes,  a result  of  0.67 
grams  is  obtained  from  the  ordinate.  Similarly  2.2  grams  of  dust  particles 
with  sizes  up  to  500  microns  have  been  collected  by  the  perfect  filter  for  the 
the  same  entry  time  of  10  minutes  and  penetration  duration  of  30  minutes. 

The  cumulative  mass  collected  for  particle  sizes  up  to  10,000  microns  repre- 
sents the  total  mass  collected  by  the  filter  since  this  represents  the  upper 
limit  of  particle  sizes  considered  in  the  analysis.  Figure  B23  presents  the 
total  mass  collected  by  the  perfect  filter  as  a function  of  time  after  weapon 
detonation  with  cloud  entry  time  as  a parameter.  Results  for  cloud  entry  times 
of  10,  18,  30,  45,  60,  90,  120,  180,  and  300  minutes  are  included  in  this  figure. 
As  an  example  of  the  utility  of  this  figure  consider  the  following:  How  many 
grams  of  dust  would  be  collected  by  a perfect  filter  in  the  aircraft  environmental 
control  system  if  it  entered  the  cloud  at  TI  = 30  minutes  and  exited  the  cloud 
after  a penetration  duration  aT  = 30  minutes?  In  this  example,  note  that  the 
aircraft  exits  the  cloud  at  TF  = TI  + aT  = 60  minutes  = 1 hour  after  weapon 
detonation.  Enter  the  abscissa  at  TF  = 1 hour  and  move  up  to  the  parametric 
curve  for  TI  = 30  minutes  and  read  from  the  ordinate  the  fact  that  1.6  grams 
of  dust  have  been  collected  by  the  filter.  The  total  mass  of  dust  collected  is 
an  extremely  strong  function  of  cloud  entry  time,  which  again  zttests  to  the 
fact  that  fallout  effects  are  significant  in  the  analysis. 

Curves  pertaining  to  the  filter  dust  mass  distribution  function,  cumulative 
mass,  and  total  mass  collected  by  the  filter  are  in  themselves  useful  from  a 
filter  design  standpoint.  However,  another  major  factor  must  be  considered. 
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The  dust  particles  collected  by  the  filter  are  a source  of  radioactivity  which 
is  carried  with  the  aircraft  and  contributes  dose  to  the  crew/electronics  in  the 
vicinity  of  the  filter.  Dose  calculations  from  this  source  have  been  made  and 
, are  based  upon  several  assumptions.  The  filter  is  assumed  to  be  a point  source 

of  radioactivity  emitting  gamma  rays  with  an  average  energy  of  1 MeV.  The 
resultant  dose  calculated  corresponds  to  a dose  “at-l-meter"  so  that  a simple 
correction  factor  of  (1  meter/d  meters)^  can  be  applied  to  the  results  to  correct 
for  actual  separation  distance  of  crew/electronics  from  the  filter. 

The  product  of  the  mass  distribution  function  and  the  specific  activity 
distribution  function  is  involved  in  the  filter  dose  calculations.  In  addition, 
after  the  aircraft  exits  from  the  cloud,  the  dust  previously  collected  continues 
to  contribute  to  the  dose  until  the  mission  is  completed.  Thus,  to  present  the 
results,  a particular  time  after  weapon  detonation  must  be  selected  as  a base- 
^ line,  and  in  this  work  the  dose  data  are  presented  at  a time  of  30  hours  after 

^ weapon  detonation.  These  results  can  be  scaled  to  other  times,  as  will  be 

demonstrated. 

Figures  B24  through  B32  present  the  resultant  perfect  filter  dose  distribu- 
tion function  as  a function  of  particle  size  and  with  penetration  duration  as  a 
parameter  for  a particular  cloud  entry  time,  TI.  The  penetration  durations 
I presented  on  each  figure  are  2,  5,  10,  15,  30,  60,  120,  and  300  minutes,  and 

each  separate  curve  presents  results  for  a particular  cloud  entry  time  in  the 
i • range  from  10  minutes  to  5 hours.  Again  the  effects  of  fallout  can  be  easily 

observed  by  a comparison  of  early  and  late  cloud  entry  data,  for  there  is  no 
^ dose  contribution  from  large  particle  sizes  for  late  cloud  entry  time.  Although 

not  as  noticeable,  it  can  be  also  observed  that  the  exponential  radioactive  decay 
I law  has  been  included  in  the  analysis,  figures  B33  through  B50. 

The  next  series  of  figures  presents  data  on  the  cumulative  filter  dose  at  30 
hours  after  weapon  detonation  as  a function  of  particle  size.  Figures  B33 
through  B41  use  semi-log  scaling  and  figures  B42  through  B50  use  log-log  scaling. 

I These  sets  of  graphs  depict  the  same  data  and  both  are  presented  to  facilitate 

' the  accurate  interpretation  of  the  numerical  values.  The  term  cumulative  filter 

dose  is  defined  as  that  portion  of  the  total  dose  attributable  to  all  particles 
collected  which  have  sizes  less  than  or  equal  to  r.  These  figures  present  this 
' cumulative  dose  for  cloud  entry  times  from  10  minutes  to  5 hours,  and  in  each 
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figure  results  are  presented  over  the  range  of  penetration  durations  of  2,  5, 

10,  15,  30,  60,  and  120  minutes.  As  an  example  of  the  use  of  these  curves 
consider  the  following  question:  For  a cloud  entry  time  of  10  minutes  and  a 

penetration  duration  of  15  minutes,  how  much  of  the  "at-1 -meter"  filter  dose 
collected  during  the  30  hours  after  weapon  detonation  is  attributable  to  parti- 
cles with  particle  sizes  of  10  microns  radius  or  less?  In  figure  B33  enter  the 
abscissa  at  10  microns,  move  up  to  the  parametric  curve  for  a 15-minute  penetra- 
tion duration,  and  read  the  filter  dose  of  6.5  Kj  rads(tissue) . The  dose 
attributable  to  all  particles  for  these  same  conditions  is  19  Kj  rads(tissue)  so 
that  the  particles  with  size  of  10  microns  or  less  contribute  approximately 
one-third  of  the  total  dose.  From  figure  841,  for  later  entry  times,  i.e.,  5 
hours,  and  for  the  same  15-minute  penetration  duration,  the  smaller  particles 
with  r <.  10  microns  contribute  1.8  K3  rads(tissue)  out  of  2.7  rads{ti ssue) . 
Figures  851  through  859,  for  the  perfect  filter  case,  provide  the  total  dose 
"at-l-meter"  due  to  sll  particles  collected  by  the  filter  as  a function  of  time 
after  weapon  detonation  and  with  penetration  duration  as  a parameter.  In  these 
figures  tesults  are  presented  for  penetration  durations  ranging  from  2 minutes 
to  120  minutes,  and  each  figure  is  for  a particular  entry  time  in  the  range  from 
10  minutes  to  5 hours.  For  the  baseline  case  of  TI  - 30  minutes  and  aT  = 30 
minutes,  note  from  figure  853  that  24  Kj  rads(tissue)  would  be  accumulated  by  a 
crew  member  located  "at-1 -meter"  from  the  filter  at  30  hours  after  weapon 
detonation  due  to  dust  collected  by  the  filter. 

This  set  of  figures,  figures  851  through  859,  for  the  perfect  filter  case 
also  provides  the  basis  for  scaling  the  previous  dose  results  (i.e.,  figures 
824  through  850)  which  were  presented  at  30  hours  after  weapon  detonation  to 
earlier  times.  A time-scaling  factor  (TSF)  can  be  obtained  from  this  series  of 
curves  and  used  to  scale  the  30-hour  results  to  the  particular  time  of  interest. 
The  time-scaling  factor  is  defined  as 


whv  .-e  D^(t)  is  the  filter  dose  at  time  t after  detonation  {tf  <.  t <.  30  hoursl, 
and  0^(30  hours)  is  the  filter  dose  at  t = 30  hours.  Both  doses  are  obtained 
from  the  appropriate  figure  (figures  851  through  859)  for  a particular  case 
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of  interest  (defined  by  a particular  TI  and  aT),  To  obtain  the  filter  dust  dose 
distribution  function  or  the  cumulative  filter  dose  at  some  time  other  than  30 
hours,  multiply  the  ordinate  of  the  appropriate  figure  by  the  time-scaling  factor. 

To  further  clarify  the  use  of  the  time-scaling  factor,  consider  the  following 
example.  Assume  that  the  aircraft  lifts  off  at  midnight  and  at  0600  it  enters  a 
radioactive  dust  cloud  resulting  from  the  baseline  massive  multiburst  dust  envi- 
ronment. Assume  further  that  the  weapons  were  detonated  at  0530.  Thus,  the 
aircraft  cloud  entry  time  is  30  minutes  after  detonation.  Assume  also  that  at 
0630  the  aircraft  exits  the  cloud.  Thus,  the  penetration  duration  aT  = 30  min- 
utes and  cloud  exit  time  is  TF  = TI  + aT  = 60  minutes  after  weapon  detonation. 
Assume  finally  that  the  aircraft  completes  its  mission  and  lands  at  a base  at 
1200  hours.  The  landing  time  corresponds  to  a time  of  t = 6.5  hours  after  weapon 
detonation.  From  figure  B53,  D.^{30  hours)  ■ ?'l  Kj  rads(tissue)  and  D.p(t)  = 

0^(6. 5 hours)  = 16  Kj  rads(tissue).  The  16  rads(tissue)  is  the  filter  dose 
(at-l-meter)  which  would  be  accumulated  by  the  crew.  The  time-scaling  function 
(TSF)  is  16  Kj/24  Kj  = 0.67.  The  6.5-hour  dose  distribution  as  a function  of 
particle  size  and  the  6.5-hour  cumulative  dose  as  a function  of  particle  size 
are  obtained  by  scaling  the  30-hour  results  in  figures  826,  B35,  and  B44  by  the 
time-scale  factor,  0.67. 

Since  crew  response  to  radiation  is  dependent  upon  the  dose  rate  as  well  as 
the  dose  received,  figures  B62  through  B70  are  included.  These  filter  dust 
rate  graphs  are  nondimensionalized  by  Kj  and  are  given  for  the  same  TI,  aT  cases 
as  the  filter  dust  dose  graphs. 

Also  of  interest  is  the  specific  activity  distribution  function,  equations 
(B-24)  and  (B-25),  which  is  shown  graphically  in  figure  B71 . 
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Figure  B12.  Perfect  Filter  Mass  Distribution  Function,  TI  = 3 Hours 


/ 


voooo 


MRS? 


MRSS  / K2  COIMeNSTONLESS) 


AFWL-TR-73-82 


Ot  Ul 

-<oi  won  u)on  uton  i/i 


PflfiTICLE  RRDIUS  IMICnaNS) 


Figure  B15.  Cumulative  Mass  as  a Function  of  Size,  TI  = 18  Minutes 


MRSS  /•  K2 


1 0000 


/ SSWW 


QOSE  / K3  : DIMENSIONLESS  DOjE  / MICRONS 


OOSE  / K3  rOlMENSrONLESS  DOSE  / MICflONJ 


DOSE  / K3  lOIMENSIOMLESS  DOSE  / MICRON 


AFWL-TR-73-82 


''I  W 

-•rv  uiC3r4  vion  ui 


PflflTICLE  molus  imOlONS) 


Figure  B30.  Dose  from  Filter  Dust  at  30  Hours,  TI  = 2 Hours 
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Figure  B31 . Dose  from  Filter  Oust  at  30  Hours,  TI  = 3 Hours 
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Figure  B32.  Dose  from  Filter  Dust  at  30  Hours,  TI  = 5 Hours 
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Figure  B48.  Cumulative  Filter  Dose  at  30  Hours,  TI  = 2 Hours 


DOSE  / K3  COiMENSIONLESSl 


AFWL-TR-73-82 


PARTICLE  flflOIUS  IMICRflNS) 

Figure  B49.  Cumulative  Filter  Dose  at  30  Hours,  TI  = 3 Hours 
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APPENDIX  C 
ZERO  FALLOUT  MODEL 

This  appendix  presents  an  investigation  of  the  limiting  case  of  a cloud 
with  no  fallout.  This  case  is  useful  because  It  yields  results  which  are  a 
worst  case  and  which  act  as  a standard  of  comparison  for  other  cloud  models 
which  include  fallout.  It  serves  as  a check  on  results  of  other  models  to 
ensure  that  the  general  trends  and  qualitative  results  are  reasonable  and 
realistic. 

This  limiting  case  is  achieved  by  “freezing"  the  10-minute  dust  density 
distribution  presented  in  appendix  A,  equation  (A-1),  for  all  time,  i.e.,  no 
later  fallout  is  allowed. 

Pd(t)  = a.  [t'^‘®  + 1.313  (C-1) 

= a^  [(0.167)’^'®  + 1.313  (0.167)‘°-^]  (C-2) 

= 19.83  a^  (C-3) 

The  same  approach  used  in  appendix  A yields  the  relationships  below.  Note 
that  the  only  time  dependence  allowed  is  that  associated  with  the  radioactive 
decay  of  the  fission  products. 

t^.  <.  t <.  tf  (C-4) 

Mj  (t) 

— 19.83  (t  - tj)  t,  <t<  t,  (C-5) 

'^1/2  ^ 


-^=  19.83 
'^1/2 


ZOO 


where  Kj,  Kj,  Kj,  and  K^  are  defined  In  appendix  A,  equation  (A-44). 

Figures  Cl  through  C7  present  representative  graphic  solutions  of  these 
dos?  rate,  dose,  and  mass  equations.  Figure  Cl  depicts  the  total  mass  of  dust 
trapped  by  the  filter  during  a penetration  of  indefinite  duration  for  TIs  of 
10  minutes  and  1 hour  after  detonation.  Note  that  the  results  are  represented 
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by  a simple  linear  function.  There  is  no  fallout  to  decrease  the  dust  pickup 
with  time.  Figures  C2  and  C3  depict  the  filter  dose  rate  as  a function  of  time. 
The  decrease  with  time  is  much  less  apparent  here  because  of  the  zero  fallout 
assumption.  The  only  decrease  is  caused  by  radioactive  decay.  Figures  C4  and 
C5  depict  the  filter  dose  as  a function  of  time.  The  entry  times  for  the  figures 
are  self  explanatory  and  the  penetration  durations  are  the  same  as  were  presented 
in  appendix  A.  Figure  C6  depicts  the  dose  rate  due  to  inmersion  in  the  cloud. 
Note  that  the  time  dependence  is  strictly  due  to  radioactive  decay.  Figure  C7 
depicts  the  cloud  immersion  dose  for  TIs  of  10  and  60  minutes  after  detonation. 
Detailed  explanations  on  the  use  of  similar  graphs  were  presented  in  appendixes 
A and  B. 

With  some  physical  reasoning  and  heuristic  argument,  this  model  could  be 
likened  to  the  cloud  resulting  from  an  atmospheric  burst.  An  atmospheric  burst 
does  not  generate  the  cloud  of  dust  as  does  a surface  detonation.  The  fission 
fragments  and  neutron  activated  weapons  debris  are  distributed  throughout  similar 
volumes,  but  there  is  little  or  no  fallout.  The  bomb  debris  consists  of  con- 
densed particles  in  the  submicron  to  micron  range  and  tends  to  remain  suspended 
in  the  atmosphere.  Therefore,  this  Zero  Fallout  or  "limit  case"  is  similar  to 
an  atmospheric  burst  and  should  yield  similar  dose  and  dose  rate  results, 
because  1n  either  ease  the  total  number  of  fission  fragments  are  the  same  and  in 
both  cases  they  remain  suspended  for  extremely  long  times  in  the  atmosphere. 

The  dose  rates  and  doses  are  functions  of  time  only  through  the  decay  of  the 
fission  products. 

For  an  atmospheric  detonation,  the  immersion  dose  determined  from  this  model 
is  directly  applicable.  The  filter  dust  mass  and  dose,  however,  cannot  be  so 
simply  determined.  In  fact,  because  of  the  very  small  particle  sizes  involved, 
very  little  of  this  radioactive  material  would  be  trapped  in  the  aircraft. 

Most  would  remain  suspended  in  the  air  and  would  be  ejected  from  the  aircraft 
with  the  air.  Therefore,  if  this  model  Is  used  to  represent  the  penetration  of 
a cloud  generated  by  an  atmospheric  detonation,  only  the  cloud  immersion  dose 
results  of  all  the  results  in  the  above  equations  would  be  pertinent. 
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APPENDIX  D 


I 

I 


COCKPIT  MASS  AND  DOSE  MATHEMATICAL  DEVELOPMENT 


The  cockpit  mathematical  development  in  this  appendix  will  assume  that 
FE{r)  = 0.  Rj,  is  the  critical  size  and  Pj(r)  is  the  probability  of  settling 
as  determined  by  the  work  in  section  V.  The  starting  point  for  this  develop- 
ment is  the  mass  rate  of  flow  to  the  cockpit  (settled  out).  In 

general , 


Mj.p(r.t)  = Mf(r,t)  P^(r) 


In  particular, 


1 < r < Rj 


R^  i r < R(t) 


Mpp(r.t) 

-V 


= 9.946  X 10- 


2 L 


1.5 


Mcp(r.t) 


= 9.946  X 10-2 


r‘0.5 


(D-1) 


(D-2) 


(D-3) 


These  equations  can  be  integrated  as  has  been  done  in  earlier  sections  and 
M(,p(r),  CMj.p(r),  and  M^p(t)  can  be  written.  Again  numerical  techniques  and 
integration  must  be  used  in  some  cases,  and  r^^  = R(tn)  ano  t^  — ''•n  — V' 

When  R(tf)  > R^' 


210 


, ,1.5 

9.946  X 10  J-r-  (tf  - ti) 


9.946  X 10'^  (tf  - t^) 


1 < r < R, 


Rj.  :£.  r <.  R(tf) 


9.946  X 10'\‘“*®  (tn  - t<)  R(tf)  < i R(ti)  (0-4) 


If  R(tf)  < R , then  these  equations  become 


.1.5 


|9.946  X 10'^  (tf  - t^) 


9.946  xlO^i  , 1.5 


1 <■  r <.  R(tf) 


<■„*  (tn'ti)  R(tf)<r^<Rc 


9.946  X 10'^  (tn  - ti)  Rj.  i S.  R(ti ) (0-5) 


Integrating  the  cockpit  mass  distribution  equations  above  over  r,  the 
cumulative  mass  equations  are  obtained.  For  R(tf)  > Rj., 


1 < r < R, 


= 3.978  x...l£l  (t  . t.)(,2.5  . T) 


(0-6) 


R(,  i !■  i R(tf ) 


= 0.1989  (tf  - ti)(r'’*® 


+ 


CM„,(Rc) 

~x~~ 


(0-7) 
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R(tf)  < r„  < R(ti) 


CMcpa(R(tf)) 

K2 

K2 

For  R{tf)  < R^, 

these  become 

1 < r <.  R{tf) 

5 1 

CP‘ 

= 3.978  X 10  - 

R(tf)  < < Rc 

CM  a(r) 

9.946  X 10"^ 

c 

R,.<r„  < R(ti) 

CM  ^(Rc) 

= — if--— + 9 

.946  X 10"^  2 ■ ti)  ir 

R(tf)  " 


(tf  - t^) 


(r"-''  - 1) 


R(tf) 


The  time  function  relations  also  yield  tv*o  sets  of  equations.  For  R(t) 

(t  - ti)^0.1989  R(tl“*®  - 3.978  x 10’^  R^,”^  - 0.1591 
R(ti) 


+ 9.946  X 10“''  r^’°'^  (t„  ' t^)  ar 


(D-8) 


(D-9 


(D-ll 


(D-r 

> R 


(D-i: 


When  R(t)  < R , 
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= 3.978  X 10"^  (t  - ti)  (R(t)^'^  - l) 


+ 9.946  X 10' 


Re'"  ft)  (tn  - t,-)  Ar 


R(ti) 

+ Z (tp  - t^)  Ar 


(D-13) 


The  cockpit  mass  equations  are  plotted  in  figures  D1  through  D19. 

The  development  of  the  dose  equations  follows  a similar  technique,  except  in 
these  equations  there  are  two  breakpoints,  and  20u,  which  R(t)  must  be  allowed 
to  progress  through.  In  addition,  provision  is  made  for  R^  to  be  less  than, 
greater  than,  or  equal  to  20u.  The  basic  development  stems  from  the  filter  oose 
rates  of  appendix  B multiplied  by  the  cockpit  efficiency  function  of  section  V, 
i.e.. 


Dcp(f.t)  = Df('-.t)  Ps(r)  (D-14) 

The  conditions  noted  above  lead  to  the  necessity  of  writing  a set  of  equations 
for  each  set  of  conditions  in  each  of  the  Dgp(n),  CD^p(r),  and  the  Djp(t) 
functions.  These  will  be  stated  below  with  little  explanation  other  than  the 
conditions  to  which  the  equations  apply.  Again  r^^  = R(tn).  where  t.j  <.  t^  <.t^. 
as  before. 

The  cockpit  dose  as  a function  of  r will  be  stated  first.  For  convenience, 
let  T(t)  = 1.25(t°-®  - t.®-®)  + ^*^i  ' ''c  - 
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If  R < 20u,  the  following  set  is  generated, 
c 


R^  < 20  < R(tf) 


R < R(tf)  < 20 


il7.64S  ^1.5 


r‘-"T(tf) 


Il7.645  T(tf) 


1 1 >■  i. 


R^  i r <.  20 


352.9  r"'-®  T(tf) 

20  < r i R(tf ) 

352.9  T(tn) 

R(tf)  < < R(ti) 

(0-18) 

17.645  ^1.5 

1 < r < Rc 

17.645  r'®"'  T(tf) 

Rj  < r i R(tf) 

17.645  T(tn) 

R(tf)  < Tj,  < 20 

352.9  T(tn) 

20  <.  Tj,  <.  R(ti) 

(D-19) 

R(tf)  < R-  < 20 


17^  ^1.5 


1 < r <.  R(tf) 


2 n 


352.9  r, 


;-®  T(t„) 

R(tf ) < r„  < Rj. 

R,  < Tp  S 20 

,-‘-"T(tn) 

20  < Tp  <.  R(ti) 

(D-20) 
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To  completely  define  i,.ie  cumulative  dose  function  the  foregoing  equations 
must  be  integrated  over  r resulting  in  several  conditional  sets  of  equations 
as  before.  When  R^.  > 2u, 


20  < Rj.  <.  R(tf) 

7,05i 
^3  = R,^ 


■T(tfj  (r"-"  - 1) 


1 < r <.  20  (0-21) 


_ CD  (20)  235. 


+ T(tf)  (r*-^  - 89.443)  20  < r < R (0-22) 


+ 705.8  T(tf)  - r'“-^) 

2^.  3..,  I 


20  < R(tf)  < R^ 


R(tf)  i.r^  iR(ti)(D-24) 


1 < r <.  20  (D-25) 


CD^„,{r)  CD^„,(20)  ,,  , ^ 

— T(tf)  (r‘*^  - 89.443)  20  < r <.  R(ti  ) (D-26) 


i!:!,  CPgJM  V r;-^T(tn)Ar 

3 h R,^  R(tf)  ^ 


R(tf)<r^<R^  (0-27) 


CD.„u('”)  C0^n,(Rc)  1 

— = —£|L_i-+  352.9  z]  T(tn)  ar  R^  < r^  i.  R(ti )(D-28) 

3 3 


I 
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iR(tf)  i20 


(r^-^-1) 


1 < r < R^  (D-37) 


+ 35.29  T(tf)  (r“‘^  - - *■  - (D'38) 


^=^^07.645  i 


CD.nu(>-)  CD,„3(20) 


S£^ 


R(tf)  <r^<Z0  (D-39) 


n 

352.9  2 T(tn)Ar  20  < < R(ti ) (D-40) 


R(tf)  ^ Rj,  <.  20 


5|!l:2..i-AaT(t,)  (r>->-i) 

^ ’’c 


S R^^  R(tf)  " 


1 < r < R(tf)  (D-41) 


CPcp3('-)  . CDcp2(Rc) 


R(tf)<r^<R(.  (D-42) 

+ 17.645  2 T(tn)Ar  R^  < < 20  (0-43) 


CD^J'-)  CPrn3(20) 


cpi*'’ '_  , rrcp3 + 352.9  2 r T(tn)Ar  20  < r.  < R(ti ) (0-44) 

^*3  MA  n (1 


20 


To  complete  the  cockpit  dose  development,  the  dose  as  a function  of  time 
must  be  determined.  This  development  again  results  in  sets  of  conditional 
equations  with  some  numerical  integration.  For  convenience,  let  TT(t)  = 

0.25  t"-®  - 1.25  t."'®  + t^.  t'“-^  . For  t^  < t < t^ , 
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20  <.  Rj  < R(t) 


(4201  .535  R ‘^  + 799. S07  - 705.8  R(t)"“-^)  TT(t) 

Kj  \ c 


R(ti) 


+ 352.9  £ r (tp  - ti)4r 

' ’ ' R(t)  " 


20  < R(t)  < R 


= R^-^  (4201.535  + 94.107  R(t)^'^)  TT(t) 


+ 352.9  - f 


RTt)  " 


(tp  - ti)4r 


R(ti) 

+ £ r (tp  - ti)Ar 

R ' 

R(t)  < 20  < Rj 


7.058  R/^  (R(t)^-®  - l)  TT(t) 

r 20 

+ 17.645  R £ r (tp  - ti)Ar 

\ 1 I G n 

_ 

Rc 

+ 352.9  R/^£  r (tp  - ti)Ar 

c 20  " 

R(ti) 

+ 352.9  £ (tp  - t^)Ar 

R,  " 


For  the  region  where  t > t,, 


AFWL-TR-73-82 


20  i R < R(tf) 


- - f°-^) 


^4201.535  + 799.907  R, 


-0.5 


R(ti) 


- 705.8  R(tf)‘°‘^)  (tf  - t])  + 352.9  £ r (t„  - ti)ir 
' R(tf) 


20  < R(tf)  < R 


(4201.535  + 94.107  R(tf ) (r^‘^)  (tf  - t,) 


R(ti) 


+ 362.9L'^E  (t„  - t^)Ar  + L (tp  - ti)J 

\ R(tf)  ^ 


D (tf) 

* K~ 


R(tf)  < 20  < R 


VII.  (t, 


7.058  R,,"^  (P-(tf)^‘^  - l)  (tf  - ti) 


20 


+ 17.G45  R, 


R(tf) 

R(ti) 

+ 352.9  S (t„  - ti)fir 

f>  '• 


t, hr.  352.9  V’§  ' 


Finally,  when  R < 20,  and  t^  <.  t t^i 
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<20  <.R(t) 

(315.634  - 28.232  R^“-^  - 7.058  R^'^  - 705.8  R(t)'“-^)  n(t) 

R(tl) 

. 352.9  E r (t„  - ti)ir  (0-51) 

' ’ ' R(t)  " 

Rj  <.  R(t)  i.  20 


= (35.29  R{t)“‘®  - 7.058  Rj‘^  - 28.232  Rj,“‘^]  TT(t) 

. (V-  - 1-»-) 

R(ti) 

+ 352.9  E r (tn  ' 

20  " 

R(t)  < Rg  < 20 


= 7.058  Rj.'^  (R(t)^-®  - l)  n(t) 

Rc 

+ (t  17.645  R E r (t„  - ti)4>- 

Vi  ' c R(t)  " 

20 

+ 17.645  E fn 
R(t^) 

+ 352.9  S r Mtn  - ti)Ar 
20  " 

When  t > t-, 


20 

117.645  E r (t„  - ti)Ar 
R(t) 


(0-52) 


(D-53) 
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The  dose  equations  are  plotted  ir,  figures  020  through  046. 
dose  rate  equations  are  plotted  in  figures  047  through  055. 

This  completes  the  mathematical  development  of  the  mass  and 
sets  for  the  cockpit  model. 


The  cockpit 
dose  equation 


CJST  MASS  I GM  / 
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Figure  D9.  Cockpit  Dust  Mass  Distribution  Function 
TI  = 5 Hours 
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FRHflCLE  flflOUB  IHICRQNS) 


Figure  D17.  Cumulative  Dust  Mass  for  Cockpit 
TI  » 3 Hours 
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Figure  D32.  Cumulative  Cockpit  Dose  at  30  Hours 
TI  = 45  Minutes 
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Figure  D53.  Cockpit  Dust  ncre  Rate, 
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